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Chapter 1 

Introduction 


In this thesis we are going to report on some of the topics investigated during the years at 
Scuola Normale Superiore, namely weak radiative decays, exemplified by the study of the so 
called Flavor Changing Neutral Current (FCNC) decays of the b quark, and weak radiative 
corrections to LEP observables, discussing the computation of the two-loop heavy top effects 
to the p parameter and to the Z —► bb decay width. 

1.1 Motivations 

These investigations contribute to a general research activity, aimed at testing the Standard 
Model at the quantum level. 

The Standard Model of electroweak interactions is a very successful theory, in the sense 
that when a process is considered which can be predicted with a good theoretical accuracy, 
there is full agreement with the experiment 1 . It appears that every known high energy 
process can be described at least qualitatively in the context of the SM; the common belief is 
that where a quantitative analysis lacks, the reason is the inadequacy of our computational 
techniques, not of the theory in itself. 

Nevertheless there are good reasons of principle to believe that the Standard Model is 
just an effective theory: for instance there is compelling evidence that p 4 models, and conse¬ 
quently the very Higgs sector of the Standard Model, is trivial. This means that these models 
cannot be renormalized unless the interaction term is put to zero: even if this result has been 
rigorously proven only for d > 4, the existing partial results speak in favour of an extension 
to d = 4. 

Leaving aside the theoretical motivations, there are aspects of the theory which are still 
essentially not verified. For instance the symmetry breaking mechanism is not clarified, and 
the Higgs sector itself lacks completely the experimental confirmation. 

These facts have motivated in the last years a considerable effort aimed at devising ex¬ 
tensions of the SM which may cure its theoretical shortcomings, without altering the “low 
energy” behavior. May be the most promising of these extensions is the Minimal Supersym¬ 
metric Standard Model. 

Every extended model predicts the existence of some new particles, but it is quite possible 
that none of these states will be detected directly in the near future: hence one has to look 


1 With some possible exceptions, for instance Rb, which is found at LEP two standard deviations higher than 
the SM expectation, a fact however that is not a serious problem at the present level of statistical accuracy. 
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for confirmations through indirect effects. The very success of the SM implies that we cannot 
expect discrepancies with the experiment at the tree level, but any SM extension may reflect 
itself in small corrections to the observables that are already being measured. 

It goes without saying that in order to disentangle new physics effects it is first necessary 
to improve the SM predictions at the level of the present experimental accuracy. In this thesis 
we shall see in two explicit examples how this program is accomplished in practice. 

The examples chosen, the b —* sy decay and the heavy top electroweak corrections, are 
good representatives of a class of studies which are expected to be particularly apt to improve 
our understanding of the SM structure: the reason is that perturbative methods appear in 
these cases well suited to extract accurate theoretical predictions. 

1.2 Plan of the Work 

The thesis is divided in two parts: in the first part we shall deal with the rare B —► X s ^ decay, 
while in the second part we shall consider the heavy top corrections to LEP parameters. 

We shall try to be self contained and to follow a pedagogical approach. The physical 
framework and the computational techniques are discussed on general grounds, as long as it 
is possible from first principles. 

In Chapter 2 we introduce from a general point of view the study of the rare B —> Ws7 
decay, which is going to be detailed in the following Chapters. We recall the recent results of 
the CLEO collaboration, which has measured both the inclusive decay and one of the exclusive 
modes, B — ► K* 7. We show how such a process is described by the Standard Model in the 
context of the Operator Product Expansion, separating the short distance contributions which 
can be computed by perturbative techniques and essentially non perturbative quantities. The 
b —► sy decay is unique in this respect, because the perturbative corrections are very large 
and enhance the decay rate by a factor ~ 4 ! This is not a sign of failure of perturbation 
theory: it only reflects the fact that the operator mixing induced by QCD corrections alters 
the pattern of GIM suppression. We shall see that despite these large corrections the inclusive 
process can be predicted accurately: terms left out in the perturbative computation are not 
expected to change the result by more than 20%. We shall briefly discuss why this process 
is interesting as a “probe for New Physics”, by comparing the predictions in the SM, in an 
extended Higgs model, and in SUSY. We shall see that one cannot expect too much from the 
comparison with experiment, unless the theoretical predictions are made more accurate. 

In Chapter 3 we start the analysis by showing how the techniques of Wilson is used in 
practice to build an Effective Hamiltonian Ti e s for rare B decays. We compute the coefficients 
of the Operator Product Expansion in the context of the Standard Model, and we start 
discussing the effect of QCD corrections on the accuracy of the matching between “full” and 
effective theory. The use of the effective theory at a mass scale much lower than the scale 
of the matching results in large perturbative logarithms, which need to be resummed (at 
least in part) to all orders in perturbation theory: this is equivalent to evolve the effective 
hamiltonian, as a renormalized operator, from the matching scale to the scale of the decay 
process in consideration. 

The techniques of “renormalization group improved” perturbation theory are briefly re¬ 
viewed in Chapter 4. We work in the context of dimensionally regularized field theories, and 
we derive the RG equation using an approach which in our opinion is more transparent than 
the usual “bare” formulation. In particular, we review some subtleties, like the so called 
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“evanescent” operators, which need to be taken into account for a consistent treatment of 
QCD corrections to weak decays. We discuss the resummation of large logarithms and we 
introduce the analysis of scheme and scale dependence: in particular we discuss how the scale 
dependence is the reflection of the error in the theoretical computation due to the left-out 
terms in the perturbative expansion. 

In Chapter 5 these tools are applied to the computation of the leading logarithmic cor¬ 
rections to b —> s'y decay. We compute the anomalous dimension matrix and use it to evolve 
the effective hamiltonian from the matching scale p — Mw to the scale of the decay p ~ mi,. 
We use a technique which has been developed to simplify the computation and which is ex¬ 
pected to be of help in evaluating the next to leading corrections. We argue why one should 
perform the next to leading computation, and how the accuracy of the existing results would 
be improved. 

In the second part of the thesis we discuss the heavy top corrections to LEP observables. 
It is known from some time that a large mass splitting in the t, b doublet results in relatively 
large O (rnf ) corrections to “oblique” and “non oblique” LEP observables, that is, related to 
self energies and vertices. Our contribution has been the evaluation of the O (mf) corrections. 
This computation can be regarded as an estimate of the theoretical error implied by the lack 
of knowledge of the full two loop electroweak corrections. 

In Chapter 6 we set up the renormalization framework, recalling on general grounds how 
it is possible to parameterize the quantum corrections to electroweak parameters. We discuss 
the basic and derived observables, and we introduce the p parameter as a way to compactly 
take into account the self energy corrections common to all the observables. We argue why 
an heavy top alters the value of the p parameter, and we also discuss the effect of an heavy 
top through vertex corrections. We briefly recall how such computations have been used to 
derive bounds on the top mass, which are in agreement with the recent result of the CDF 
collaboration (mt — 170GeV). 

In Chapter 7 we show how the heavy top corrections can be evaluated by using a reduced 
model, the so called “gaugeless” limit of the Standard Model. Indeed, in the limit of gauge 
couplings much smaller than Yukawa couplings, we can forget the propagation of gauge fields, 
and the model needs only to include the “would be” Goldstone bosons, which are coupled to 
fermions through the Yukawa term. We show how this limit is well defined starting from the 
Standard Model formulated in the background field gauge, and we show how the self energy 
and vertex corrections needed to evaluate the p parameter and the Z —> bb decay width are 
related by Ward Identities to correlation functions evaluated in the reduced “Yukawa” model. 
We show how this Ward Identities reflect the conservation, in the reduced model, of the 
currents coupled to the Z and W ± fields. 

In Chapter 8 the computation of the heavy top effects is detailed: we define the renormal¬ 
ized model by introducing the subtractions needed to make all the results finite, and we clarify 
the connection with the basic observables. We compute the one and two loop corrections and 
we discuss their numerical significance. Then we consider the effect of the corrections on 
some of the LEP observables. We discuss the validity of the approximation and we argue 
why one should not consider the heavy top effects as the major source of correction to LEP 
observables. 

The conclusions (Chapter 9) will be devoted to resume the results obtained and to discuss 
further possibilities, in particular the computation of the Next to Leading Corrections to B 
decays we are currently involved with. 

In the Appendixes we have collected some useful formulas. In particular, in App. A we 
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fix the notations, relating pseudo-euclidean and minkowskian metrics. In App. B we give 
definitions and identities for the Dirac algebra in d dimensional space. In App. C we briefly 
recall the Standard Model formulation, its quantization and the relevant Feynman rules. 
Finally in App. D we briefly recall the background field gauge method and its application to 
the Standard Model. 




Part I 

The B —» Xg'y decay 



Chapter 2 


Overview 


2.1 The CLEO results 

The CLEO collaboration, at the Cornell Electron Storage Ring (CESR), has reported the first 
observation of rare B decays connected at the quark level to the b —> s'y process [ 88 , 105]. 

At Cornell B, B mesons are produced at the T — 4S 1 resonance, a J = 1 unstable state of 
mass m = 10.5800 ± 0.0035GeV and total width T = 23.8 ± 2.2MeV. Both neutral B°, B° 
and charged B~, B + pairs are produced, accompainied by light quark pairs, uu, dd, ss, cc, 
whose background can be rejected by subtracting measures taken about 50MeV below the 
resonance. 

The B mesons are produced with a cross section ratio to the continuum roughly of 1 to 2.5, 
and decay mainly in charmed mesons, through Charged Currents interactions. These modes 
are used to determine the mass of the B meson and the results show the good resolution 


M b o = 5280.3 ± 0.2 ± 2.0 

Af B - = 5279.9 ± 0.2 ± 2.0 . (2.1) 

The CLEO collaboration has collected [105] 2.15 million B meson pairs at the resonance, 
and about 6.6 million below resonance. Already in a smaller sample, ( 1.2 million pairs) 
the collaboration has reported [ 88 ] the observation of 13 events corresponding to the rare 
B° A"*(892)°7, B~ —> A"*(892)“7, resulting in a branching fraction 

BR (B -> AT* 7 ) = (4.5 ± .5 ± 0.9) x 1(T 5 . (2.2) 

These exclusive decay modes were already expected to be dominant among the channels 
relevant for the b — ► sj transitions, since the A '7 decay was excluded by angular momentum 
conservation, and the K*( 892) is the lightest vector meson resonance allowed in the final 
state. 

The two body nature of the decay permits a rather clean reconstruction: the experimental 
signature is an hard photon (2.6GeV) recoiling against the decay products of the excited A* 
state, a K ir pair. The K n pair mass is reconstructed and required in a range 821 < Mk n < 
971MeV, and various cuts are imposed to reject continuum [ 88 ]. 

It has also been possible to measure the inclusive decay rate, BR(B —> X s ~f), by looking 
for an hard photon without reconstructing the accompaining hadronic state but requiring its 
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Figure 2.1: One loop diagrams for b —» sj in the SM 

mass in the typical strange meson region, 0.5 — 2GeV. The signature is a photon with energy 
between 2.2 and 2.7 GeV, this spread resulting from the Fermi motion of the b quark inside 
B meson. Referring to [105] for a discussion of the identification and background rejection 
techniques, we can quote the experimental result 

BR (b -> 57 ) = (2.32 ± 0.57 ± 0.35) x 10 ~ 4 ; (2.3) 

we shall see that this observation is quite consistent with the Standard Model predictions. 
It is also useful to observe that the following inequalities are true at 95% confidence level 

1.0 x 10 “ 4 < BR (b —> ai) < 4.2 x 10~ 4 , (2.4) 

a range of variation important in order to constrain all the non-standard model contribu¬ 
tions to the decay. 

To size the importance of this rare decays and of the hard experimental effort needed to 
detect them, let us consider their theoretical interpretation. 

2.2 The structure of the B —> X <,7 decays 

The b —> sy decay, a Flavor Changing Neutral Current process, is not generated in the Stan¬ 
dard Model classically, that is, at the tree level. This feature is “protected” at the quantum 
level by the GIM mechanism: in the limit of equal masses for the up-quarks, the diagrams in 
Fig. 2.1 cancel exactly, and the flavor diagonal structure of the Neutral Current interactions 
is preserved from receiving quantum corrections. 

As soon as it was recognized that the t quark mass had to be rather large [ 11 ], it became 
apparent the violation of the GIM mechanism, and the generation of FCNC effective inter¬ 
actions, which can be accounted for in an effective low energy theory which incorporates the 
contributions resulting by the virtual effect of the heavy top. 

We shall discuss in the next chapter how the effective theory is generated, for the time 
being it suffices to observe that thanks to the large mass of the W boson, in analogy to the 
Fermi theory, decays like B —> K *7 can be described at short distances, by a local magnetic 
momentum operator 

g 

Cm.m. (m-ft) — “2 , 


(2.5) 
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with a coefficient C, which depends on the scale of the process and on the short distance 
physics details. This parameterization is specific to the SM only in the fact that the unitarity 
of the CKM has been assumed, and the V ub V* s have been dropped, but in any effective 
low-energy theory the essential characteristics do not change. 

Note the explicit nib factor: in the limit of mb —► 0 the b helicity would be conserved and 
there would be no way to emit a real vector particle. The amplitude for B —► Ji *7 is given 
by the on-shell matrix element [82] 

A = (m b ) V tb VZrf{K* \sa^b R \ B) ( 2 . 6 ) 

where q, r 7 are momentum and polarization of the outgoing photon. The matrix element 
contains the long distance physics, the details of the hadronization, and can be parameterized 
as follows [64] 

(K* (k, e) \sa l _ w q u b R \ B(p)} = (q 2 ^j + C 2 T 2 (V) + terms zero on-shell (2.7) 

where q = p — k and 


CA = e fl „x p s u p > 'k p 

Cl = (m 2 B - m**) - £ ■ q (p + k)^ . (2.8) 

On shell T 2 ( q 2 = 0) oc T\ (q 2 = 0), and after proper phase space integration and sum over 
polarizations one gets 


r {B -»• K* 7 ) = 


a 


32t r 4 


m b G 2 F m% 


1 - 


m K* 


m: 


\V tb V*\ 2 \C(m b )\ 2 


B 


T 


q 2 = 0 


(2.9) 


This expression deserves a few comments: 

• the rate is proportional to oqedG^, while most other FCNC processes involving leptons 
(like b —> X s e + e~) or photons are of order Uqed C 2 p- 


• The CKM coefficients appearing in the expression have not been directly measured, but 
in the SM they can be deduced from unitarity, and are known to be 


| Vtfo| = 0.9987 -0.9994 

\V ts \^\V cb \ = 0.044 ± 0.006. (2.10) 

If one considers a related FCNC process, the B —► X ^7 decay, one immediately deduces 
that 


BR{B^X dl ) \V td \ 2 
BR{B-*X sl ) |I 7 S | 2 


( 2 . 11 ) 


which is a consistent suppression. 
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• The separation between short and long distance physics is made, in the context of 
OPE, at a scale m b ~ 5GeV which is well within the perturbative QCD range. We shall 
elaborate on this subject in the following chapters, but the essential is that the large m b 
mass makes us confident that the factorization of short and long-distance contributions 
is reliable. 

• The long-distance physics, contained in the Tj form factor, has been considered for a 
long time the larger source of uncertainty. 

Recently however, as pointed out in [104], a remarkable agreement between Lattice QCD 
results and QCD sum rules has shown up. 

On the lattice two independent groups have computed the form factor T\ (q 2 = 0); Bernard 
et al. [64] obtaining j' 1 Bernard = 0.20±0.02±0.01, Bowler et al.(UKQCD collab.) [82] obtaining 
TjUKQGD = 0.30+}, 0 . 

From QCD sum rules other results are available: j , i CoIdllgcl ° ct aL = o.35 ± 0.05, obtained 
in [84], 7j Ba11 = 0.37T0.05 [83], T 1 Narison = 0.31±0.013±0.06 [101], T^ li et aL = 0.32T0.05 [97], 

This results determine the ratio 


R k * 


BR (B -> K* (892) 7 ) 
BR (B - X sl ) 



( 2 . 12 ) 


of the esclusive over inclusive process, and while the lower Lattice results gives a value of 
the order of 6 %, the larger Lattice and the QCD sum rules results give a value larger than 
20 %: despite the large errors, this result confronts well with the experimental result 


R \= 0.19 ±0.09 . 


(2.13) 


The agreement is quite encouraging, and calls for further study in order to improve the 
confidence in the long-distance physics analyses, and eventually reduce the difference in the 
numerical results. 


2.2.1 The inclusive rate 

The inclusive rate F (B X s 7 ) can be accurately predicted on the basis of perturbative QCD 
only. The reason is that in a spectator model the probability of the decay is determined solely 
on the basis of the amplitude for the quark process b —> sy, that is, one assumes that no 
interferences are present between short and long distance effects. 

This is more than an assumption: the Heavy Quark Effective Theory allows to demonstrate 
that the spectator model is a well definite limit of QCD, the leading term in an expansion 
in powers of - 7^7 (see for instance [100] and references therein). Moreover it can be shown 
that the expansion starts at the A 2 / m Q level, hence for the b quark, mb — 5GeV, the non 
spectator corrections are expected to amount to less than 1 %. 

In the spectator model the inclusive rate is therefore given by 

L (B X s7 ) = a ^J m b \V tb V t *\ 2 |C m . m . (m b )| 2 ; (2.14) 

a more refined discussion of this formula, and of the dependence on the scale used to 
separate short and long distance contributions shall be given in the following chapters. 
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The predictivity of this formula is enhanced by reducing the strong dependence on the 
quark mass mb, normalizing the rate to the semileptonic b decay, T (b —> ceV e ). The ratio 
results 


T (6 —■> s 7) 

T (6 —* ce P e ) 


\VtbV t *\ 2 6a Q ED ! 

^ '('-/(Sf))' 


(™ b )| 2 , 


(2.15) 


where the / function results from the integration over the phase space in the semileptonic 
b decay 


/ (x) = x 2 ^8 — 8x 4 ± x 6 ± x 2 lnx) . 


(2.16) 


The ratio in Eq. 2.15, or equivalently the ratio of branching fractions 


BR [B -> X sl ) 
BR{B^X cl ) 


(2.17) 


together with the experimental measure [70] 


T(B -> X c eV e ) = 10.7 ±0.5 


(2.18) 


allows to predict the branching fraction, slowly dependent on mt, shown in Fig. 2.2, which 
includes the short distance QCD corrections at the level of Leading Logs. The error bars 
accompaining the theoretical values have been computed following [81], and more details will 
be given in Sec. 5.4.1. 

Despite the large errors, the figure shows a remarkable success of the SM and of pertur¬ 
bative QCD. We shall see indeed that the theoretical prediction is strongly influenced by the 
aforementioned short distance QCD “corrections”, and it is enhanced by a factor ranging 
between 4 and 5 for reasonable values of the t mass. It is worth recalling that even before 
the measurement of the inclusive decay rate, the knowledge of the BR{B —> K* 7 ) and the 
SM prediction for the inclusive rate (including the short distance QCD corrections) allowed 
to obtain a value for the ratio defined in Eq. 2.12 


R k * = 0.15 ±0.06 . (2.19) 

In absence of the QCD enhancement this ratio would be of the order of 60%, showing 
how dramatic is the effect of these corrections. One may wonder if so large corrections can 
be reliable, and in which sense they are “perturbative”. 

2.2.2 The QCD enhancement 

The substantial short distance QCD enhancement makes the b —> s'y process the only known 
SM process “dominated by two loop effects” [81]. 

Let us briefly comment on the reasons of this big enhancement. 

As already stated, at short distances the Feynman diagrams in Fig. 2.1 generate the 
decay. This graphs give rise to a coefficient C m m for the magnetic momentum operator in 
Eq. 2.5, whose value varies between ~ 0.17 — 0.21 for a top mass in the interval 140 — 200GeV. 
However, through QCD corrections also the Fermi type operators, like ( sc) v _ A ® ( cb ) v _ A , 
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which appear in the Effective Theory below the Electroweak scale, contribute to the process 
via two-loop Feynman graphs, as in Fig. 2.3. 

These operators give rise, order by order in the loop expansion, to contributions propor¬ 
tional to powers of as (Mw) log — 0.33, where the value of the fj, scale has been fixed by 

matching the effective theory and the underlying complete theory (say, the SM). This correc¬ 
tions can be summed in part to all orders of the perturbative expansion, and can be included 
in a redefinition of the coefficient of the magnetic momentum operator 1 (schematically) 

( cx [Tin) \ {ex (m) \ 

a s \M w )) Cmm - (Mw) + -" + fF ( a s \ Mw )) CVermi ( 2 ‘ 20) 

where the m argument of the coefficient means that the renormalization scale appropriate 
to compute the matrix element in Eq. 2.6 has been “scaled down” to m, and the / functions 
contain the QCD corrections, which can be expressed in terms of the ratio of the QCD 
couplings at different scales. 

We shall see that the additive renormalization schematically addressed by the “Fermi” 
coefficient in Eq. 2.20 is responsible for the consistent improvement, while the multiplicative 
renormalization tends to suppress the rate. 

The physical origin can be traced back to the evasion of the GIM mechanism. If the t and 
c masses where equal, the diagrams in Fig. 2.3 would have been canceled by corresponding 
graphs with a virtual t quark. 

2.3 The b —> sy as a probe of new physics 

To be definite, in Sec 2.2 we have briefly addressed the origin of this rare decay in the context 
of SM. A more complete discussion will be given in the following chapter, but it should 
be stressed that the b —> s'y process is not only a test of the SM but more importantly is 
a window in the higher energy structure of the theory. The reason is that non standard 
contributions resulting from the exchange of undiscovered particles are comparable in size to 
the SM contributions, which are suppressed by the symmetry mechanism preventing large 
FCNC effects. Hence on general grounds one expects that even if it will be not possible for 
some time to extend the investigation of the spectrum of elementary particles, the indirect 
effects will be accessible through the combined study of the different rare decays. 

Fixing the attention to b —> 57 , many detailed analyses exist, in the context of the Minimal 
Supersymmetric Standard Model (MSSM) [59, 63, 86 , 99], the Two Higgs Doublet extension 
of the scalar sector [79, 80, 81] and in other models, as in the SU(2)l x SU(2)r x 17(1) 
extension of the electroweak gauge group [89]. 

The lack of space and direct experience prevents us from adequately reviewing the many 
interesting contributions to the argument, but a few general considerations are possible. 

2.3.1 Two Higgs doublet models 

The two Higgs doublet models can be seen as the mildest extension of the electroweak sim- 
metry breaking sector of the SM, and the resulting charged Higgs fields in the spectrum 
naturally induce FCNC transitions. There are two versions of the model, the so called Type 
I in which one of the doublets “gives mass” to all the fermions while the other decouples, and 


1 In the leading log approximation 
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Figure 2.4: Charged Higgs contribution to b —> 57 . 


the Type II model in which one Higgs doublet (^ 2 ) couples to up-quarks, while the other (yy) 
to down-quarks. The Type II model is more popular because it may possibly give a “natural” 
explanation of the large mass splitting in the t, b doublets, in terms of a large ratio of the 
vacuum expectation values, commonly parameterized as 

Vo 

tan/3 = — ; (2.21) 

Vl 

moreover, this model appears as the minimal extension of the Higgs sector when consid¬ 
ering SUSY theories. 

In either models there are new contributions to b —> sy through the diagrams in Fig. 2.4, 
with the general form [80] 


Hsm 


2 

m t 
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( 2 . 22 ) 
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where Hsm and A£ are the SM and charged-Higgs contributions respectively; in model 
I one has A = — 1/tan 2 /3, while in model II A = 1. For small values of tan/3 both the models 
give an enhancement of the amplitude, and in model II the term A\j gives an amplitude 
always larger than in the SM. However in model II large tan /3 values are more appealing, as 
discussed above. 

In Fig. 2.5 we show the prediction for two values of Mh±, already considered in [81]. 
For a large Mh ± the SM and the 2HDM become undistinguishable, at the present level of 
theoretical accuracy. 


2.3.2 b —> sy in the MSSM 

Good reviews on the subject exist, and we shall limit to a brief account, mainly following [59, 
78, 99]. Even restricting the attention to the so called N = 1 models, that is, the SUSY 
models which are derived as low energy theories from a N = 1 supergravity model, a large 
number of parameters exists. The MSSM is obtained with some additional constraints 

• a new symmetry, the R parity, is imposed in order to avoid large baryon-lepton number 
violations. 

• Only the minimal number of fields needed to “supersymmetrize” the known spectrum 
are introduced. 
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Figure 2.5: B —■> X s 7 in the 2HDM 

• Canonical kinetics terms are chosen for all the scalar fields. 

While the particle content is fixed, with these choices there are essentially 4 new parameters 
with respect to SM (5 if one relaxes the assumption of a flat Kalher metric which underlies the 
choice of canonical kinetic terms), which determine the low-energy theory. By requiring also 
that the correct scale of electroweak symmetry breaking is reproduced, and not introducing 
a third Higgs doublet, the parameters reduce to 3; of these, one is again commonly chosen as 
tan f3, while the other two can be for instance chosen as two mass scales characterizing the 
pattern of SUSY breaking. 

The enlargement of the spectrum, with respect to the SM, brings naturally new sources 
of FCNC contributions: apart the extension of the Higgs sector, discussed in the preceding 
section, one has an immediate new contribution coming from the supersymmetrization of the 
W and charged Higgs contributions. The fermion mass eigenstates resulting from W and 
H± mixing, the so called charginos x±, are exchanged in loops where the up squarks also 
circulate. 

Other peculiar contributions result from the Flavor Changing vertices connecting quarks 
and squarks through the exchange of neutralinos x° and gluinos. In summary one has con¬ 
tributions resulting from the following virtual particles 

1. W~ and up quarks, 

2. charged Higgs fields H~ and up quarks, 

3. charginos, x~, and up squarks, 









2.3 The 6 —> sy as a probe of new physics 


19 


4. neutralinos, x° and down squarks, 

5. gluinos, g , and down squarks 

where the last three are peculiar of SUSY. In some of the older analyses the first two 
contributions were considered dominant, thus giving bounds analogous to the one considered 
in the SM. However more complete studies [59, 86] show that only the neutralino and gluino 
contributions can be numerically neglected. On the other hand, while as in the 2HDM the 
charged Higgs contributions interfere constructively with the SM contribution, there exist 
large portions of the parameter space in which the chargino amplitude interferes destructively, 
and it can also become the dominant contribution for a large tan (3 value. 

It is not possible even to account for the exploration of the parameter space needed to 
make these statements quantitative, so we refer to a recent analysis [99]. 

We can expect in a near future an improvement of the accuracy of the experimental results, 
however we see that the possible interference effect among different SUSY contributions makes 
difficult a clear distinction from the SM. This is even more true in presence of the rather large 
theoretical uncertainties, mainly originated by the lack of knowledge of higher order QCD 
corrections. 




Chapter 3 


Effective theories 


3.1 Generalities 

In this thesis we consider physical processes which are strongly influenced by quantum cor¬ 
rections. A direct consequence, relevant both for weak decay processes, and the electroweak 
corrections to LEP observables, is the dependence of the theoretical predictions on widely 
different mass scales. 

The reason is that the SM 1 has the ambition to account, on quantitative grounds, for 
phenomena occurring at largely different energies, starting from a common fundamental de¬ 
scription. For instance, in many weak decays the typical energy scale, the Q value of the 
reaction, is of the order of lGeV, while the fundamental interactions result from the exchange 
of virtual particles with masses of the order of lOOGeV. On the other hand, observables 
determined with experiments at energies of the order of the Z resonance, must be related to 
low-energy experiments, like the /3 decay or the v N scattering. 

A better understanding of the phenomena, as well as a practical way to properly ac¬ 
count for the different scales involved, can be obtained by using the effective field theory 
formalism [26]. 

3.1.1 What is an Effective Theory? 

A number of useful reviews exist (see for instance [27]), so we shall limit ourselves to a very 
general discussion, and then proceed through the study of a practical case. 

The idea is quite simple: we start from an high energy theory, say the SM, described by 
a local lagrangian 2 . Suppose it involves light fields, cj), and heavy fields, T. 

We are going to study a process which can be described in terms of correlations of light 
fields 

(0 \<fi(xi)c/>(x 2 ) ■ ■ ■ 4>(x n )\ 0) = f dtdte-SoWKxMxi) ■ ■ ■ <j>{x n ) . (3.1) 

We can perform the functional integration over the heavy field, and have an equally good 
description of the light fields interactions. The same Green function will be obtained from a 

x Or any sensible extension of it. 

2 In some sense, this is the definition of a fundamental theory, at least for the majority of physicists. 
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functional integral 

J d^e~ Sl ^^{xi)^{x 2 )...(t){x n ) (3.2) 

involving an effective action Si, which is not, at least in general, local. If we try to express 
the action in terms of local interactions, we find that an infinite tower of composite operators, 
involving more and more elementary fields, is required in order to match the two descriptions. 

Technically this procedure, which allows to go from an effective non local action to a local 
theory is the Wilson’s Operator Product Expansion: in presence of an heavy mass-scale, it 
allows to organize the infinite tower of non-renornralizable interactions as a Taylor series in 
inverse powers of the heavy field mass. The short distance physics is then “buried” in the 
coefficients of the resulting local operators, while the long distance physics remains explicit. 

It is worth noting the similarity with the description of critical phenomena and the exis¬ 
tence of universality classes: long distance physics is described by an entire class of different 
theories, differing on the amount of short distance physics which is left explicit. 

3.1.2 Matching 

In all the practical cases the non-local effective action cannot be determined: instead, one 
builds a local effective lagrangian, or more generally an effective Hamiltonian, up to some 
power p of the inverse mass (for instance some order in Gf ), which will involve a finite 
number of operators and coefficients. 

For instance, in the case of weak decays, for any initial and final state, which can be 
written in terms of “light” fields, it must be imposed the indentity 

</K)f ull M = </ l^effl i) (/*) + O (3.3) 

where H e s is a sum of local renormalized operators N [O], 


H eS = J2Cj(M, p)N[Oj\ (3.4) 

j 

built of light fields, and depending on M through the coefficients C. The bracket on the 
left hand side addresses an amplitude i —> f evaluated in the “full” theory, while on the 
right hand side the expectation value of the effective hamiltonian between the same states is 
evaluated, as specified by the lagrangian of the light theory 3 . 

A matching scale p is therefore introduced, as the scale where coefficients of the low-energy 
theory are tuned to reproduce the results of the high energy theory, on the basis of a finite 
number of test processes. 

On physical grounds, the matching scale can be regarded as separating short distance from 
long distance physics. In fact, through the expansion in loops both theories generate loga¬ 
rithmic corrections, depending on the typical momenta and masses involved and an arbitrary 
renormalization scale. This means, as well known from QCD, that choosing the renormaliza¬ 
tion scale of the same order as the typical scale of the process, the higher order corrections 
are minimized. 

3 A1so the parameters of the light theory undergo finite renormalizations in the matching procedure, but we 
shall see that they are unobservable, and therefore irrelevant. 
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In the matching procedure, one is guided by the “complete” theory to set the fi scale in 
order to have theoretical predictions for the test processes as accurate as possible: so the 
effective theory is determined at a large scale, while it is to be applied at processes at a lower 
scale. 

3.1.3 Renormalization Group 

The coefficients of the effective theory, determined at the matching scale, say n — M , act as 
the experiments setting the parameters of a renormalizable field theory. As in a renormalizable 
field theory, one determines the fundamental parameters, the “couplings” of the theory, at 
some scale, and then the theory is used at a possibly largely different scale. 

Now, at a fixed order in the expansion the effective theory is renormalizable, and the 
coefficients can be scaled using the renormalization group evolution, in order to improve the 
convergence of the perturbative expansion. 

It is worth then reporting the observation of Georgi [48], “the renormalization group is 
simply the matching of the theory at the scale fi to the theory at the scale n — dfj,, without 
changing the particle content”. In other words one integrates out the effect of the energy 
modes higher than the process in consideration. 

The RG evolution results then from the integration of this differential matching procedure, 
and allows to use in a optimal way the informations coming from perturbation theory. 

The Effective Hamiltonian is scaled down from the matching scale // ~ M to the charac¬ 
teristic low energy scale of the physical amplitude 

- M ’ M ) -> Y^C(n ~ m, M) {f\N[Oj}\i) (/x~ m)(3.5) 

3 3 

which means that one determines the coefficients appropriate to the operators renormal¬ 
ized at the scale // ~ m. 

In this way, the large logarithms that would appear in the matrix elements as a conse¬ 
quence of the QCD corrections in presence of two widely different scales, the renormalization 
scale of the operators and the typical scale m of the low energy process, are transferred in 
the coefficients. 

In the rest of the Chapter we shall give more details on this procedure, illustrating the 
matching procedure for the Effective Hamiltonian which describes the rare b decays. 

3.2 The FCNC Effective Hamiltonian 

In the Standard Model, our “fundamental” theory, the rare b decays are mediated by the 
graphs in Fig. 3.1, where in the loop circulate the c and t quarks 4 , and the virtual dashed 
lines are the W boson or the charged component of the Higgs multiplet. 

A particular gauge is chosen, the so called R £ gauge, described in Appendix C: in this 
gauge the theory is renormalizable by power-counting, and the nonlinear form of the gauge 
fixing term cancels trilinear vertices A, W, q i, reducing the number of graphs to be computed. 
Most importantly, thanks to the covariant derivative in the gauge fixing term, the electromag¬ 
netic Ward-Takahashi identities are unbroken, and there is no need of finite renormalizations 
to reinforce them. 

4 the u quarks can be completely neglected thanks to the smallness of the corresponding CKM matrix entries 
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Figure 3.1: b —* s'y process in SM 



Figure 3.2: Contribution to Effective Hamiltonian from the t quark 


The first step in the determination of the effective hamiltonian is the computation of 
the amplitude for a minimal set of relevant processes in the full theory. In the following we 
will integrate out the t and W, ((f>) fields: graphs with heavy particles t, W give rise, in the 
expansion in powers of external momenta, to local interactions which can be reexpressed as 
a sum over local operators. This is pictorially expressed in Fig. 3.2 

By computing the proper vertex part, one finds easily the following effective hamiltonian, 
resulting from the t exchange only 5 : 

H eff :top = ^- V * s V tb J 2 C \ O t (3.6) 

v 2 1=1 


where 


5 Recall that = f -^=j jir- is the coupling of two weak currents 
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0\ = 0 SL £)DDb L 

(47 T)- 

O 2 (47t)^ ^ ^ Qd)SL ^ F - CT^jy ^ ft/, 

O3 = (4yj-^2 ( ^ 6 Qd)SL TV bljD 

O4 = - 1 » rrib s L DD b R 

(4tt) z 

0 5 = e Qd) m b SL b R 

° 6 = (4^)2 (*&>) SL {A T A Gf lv o>} &L 

Or = j^{-ig s )s L luT A b L {D lx G lxv ) A 

° 8 = (4^)2 (- i 9s) m bSLT A G A I/ (T IJll/ b R (3.7) 

and the computation gives 


A 

A 

A 

A 

r l 

°5 

A 

A 

A 


(2 + x) (1 — 5x — 2x 2 ) | x 2 (2 + x) log(x) 

6 (—1 + x) 3 + (-1 + x) 4 

46 — 141 x + 105 x 2 + 8 x 3 (2 — 5 x) x (—2 + 3 x) log(x) 

24 (-1 + x) 3 4 (—1 + x) 4 

104 — 312 x + 237 x 2 — 47 x 3 (8 — 32 x + 54 x 2 — 30 x 3 + 3 x 4 ) log(x) 

36 (-1 + x) 3 6 (-1 + x) 4 

x (l — x 2 + 2x log(x)) 

2(1 - x) 3 

(3 — 5 x) x x (—2 + 3 x) log(x) 

4 (-1 + x) 2 2 (-1+x) 3 

—8 + 12 x + 15 x 2 — x 3 (2 — 5 x) x log(x) 

24 (—1 + x) 3 4 (-1 + x) 4 

—4 — 42 x + 21 x 2 + 7 x 3 (—4 + 16 x — 9 x 2 ) log(x) 

36 (-1 + x) 3 6 (-1 + x) 4 

(—3 + x) x x log(x) 

4 (—1 + x) 2 + 2 (—1 + x) 3 


with x = mA 

This is not the end of the story, because the graphs with charm exchange cannot be 
accounted for in the same manner: in fact, one easily recognizes that coefficients C 3 and C 7 
diverge for x —► 0 , hence the limit m c —► 0 in the full theory leaves graphs which are infrared 


divergent in the limit of zero external momenta and cannot be expanded in series. This is 
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Figure 3.3: Contribution to Effective Hamiltonian from the c quark 


what we expect, because of the contribution, in the light theory, of the Fermi-type interaction 
terms, resulting as the low energy limit of the tree level W-exchange graphs: this interaction 
is to be included in the effective hamiltonian as 


4 G f 


VtsVtbOl 0 = 


^cs^cb s oTlfj, c /3 ® ba 


s -,Li 




(3.9) 


where the unitarity of the CKM matrix V has been used; a, (3 are the color indices, and 
the 9 suffix has been reserved for an operator which will be generated by QCD corrections, 
having a Fierzed color structure 

s a ^c a ®c^b 0 . (3.10) 

The 4—fermion operator gives rise to the infamous “penguin” graph, in the right hand 
side in Fig. 3.3: the difference between the graphs in the full theory, and this “penguin” graph 
is local and can be easily computed. It does not contain IR divergences, because the full and 
effective theory have the same long distance behaviour, hence it can be expanded in powers 
of the external momenta. It results the “charm” contribution to the coefficients in Eq. 3.8, 
again taking into account the unitarity of CKM matrix 
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(3.11) 


it is important to note that the coefficients C 3 , C 7 are explicitly dependent on the subtrac¬ 
tion scale Hi as we expected from the general considerations: this dependence is completely 
canceled by the corresponding dependence in the matrix element of the operator 0\q. This 
is true, at this level, by construction, but we shall see that the inclusion of QCD corrections 
shall make the cancellation only approximate. Moreover, they appear pole parts in the co¬ 
efficients: the reason is that the matrix element of operator 0 10 is divergent, while the full 
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theory amplitude was finite. This divergence is canceled by the coefficients, hence we can 
redefine the effective hamiltonian as made of renormalized operators, for instance in the MS 
scheme. 

Therefore we have the effective hamiltonian 

H eS = -^=-VtsV tb ^2 C*N [Oi\ (3-12) 

v 2 i=1 

with the corrected coefficients 

x (1 + 5 x) x 2 (2 + x) log(x) 

2(1 - xf + (-1 + x ) 4 

x (—1 — 11a; + 18x 2 ) (2 — 5x) x (—2 + 3x) log (a;) 

8 (—1 + x ) 3 4 (—1 + x ) 4 ’ 

—16 + 48 x — 73 x 2 + 35 x 3 (8 — 32 x + 54 x 2 — 30 x 3 + 3 x 4 ) log(x) 

12 (—1 + x ) 3 6 (—1 + x ) 4 

x (1 + x) x 2 log(x) 

2 (—1 + x ) 2 + (1 - x ) 3 ’ 

(3 — 5x) x x (—2 + 3x) log(x) 

4 (—1 + x ) 2 2 (—1 + x ) 3 

(4 — x) x (—1 + 3x) (2 — 5x) x log(x) 

8 (—1 + x ) 3 + 4 (—1 + x ) 4 ’ 

8 — 42 x + 35 x 2 — 7 x 3 (—4 + 16 x — 9 x 2 ) log(x) 

12 (—1 + x ) 3 6 (—1 + x ) 4 

(—3 + x) x x log(x) 

4 (—1 + x ) 2 + 2 (—1 + x ) 3 ’ 

0 

1 . (3.13) 

We stress that the results quoted are obtained with operators renormalized in the MS 
scheme, and the NDR scheme for the treatment of 75 : when considering QCD corrections we 
shall see how to compare different schemes and to ensure the scheme independence of physical 
results. 


Ci = 

C 2 = 
C 3 = 

C'A = 

C 5 = 
C 6 = 
C 7 = 

C 8 = 

C 9 = 
C10 = 




Chapter 4 

Renormalization Group improved 
perturbation theory 


In this somewhat technical Chapter we shall give a brief overview of perturbative renormal¬ 
ization, avoiding the language of “bare” fields and operators; we shall rederive some of the 
formulas of Renormalization Group, which shall be used in the following Chapters. 

4.1 Notations for the Dimensional Regularization 

Let us consider a field theory regularized by continuation in the number of dimensions {d = 
4 — 2e), and described in lagrangian formulation by the action 

S = Jd d z (V y £ (. 9i + Pi) Pi (*) (4-1) 

i 

where the lagrangian density is intended as a sum over local operators £j, each multiplied 
by a charge gt , and Pj denotes the counterterms needed to render the Green functions finite 
in the e e-> 0 limit, and also to impose particular renormalization conditions. Note that we 
suppose that every operator, including the terms in the quadratic part of the lagrangian, 
carries an explicit coupling, which will be very useful to deduce general identities. We shall 
also use the abbreviated form 

Si = j d d z (/)~ £ U . 

Further, let us consider a composite operator Gj = E^Oi , where Oj is a product of some 
of the fields in the lagrangian, and Ei is a product of couplings. The renormalized product is 
then defined as 

N [« Gt} = N [EjOi] = EiMijOj , (4.2) 

j 

where the sum is carried over a number of operators (a complete basis) sufficient to 
guarantee that with an appropriate choice of the ” counterterms” matrix M every Green 
function with the insertion of the operator Gi is made finite in the enO limit. 
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In this formulation, the mass parameter appears in the lagrangian explicitly to compensate 
the change in mass dimensions of the integration measure, and therefore fields and couplings 
have their canonical (four dimensional) mass dimensions. 

We now want to determine the change in the couplings in order to compensate for a change 
in the scale g. 


4.1.1 The renormalization scale dependence 

Let us consider for definiteness a Green function, the correlation of fields and of a 

local operator O , and take the derivative with respect to /j, 


/r^(O|iV[G i ]0i...0 n |O) = /o 


N[G l ]2eJ2(9i + Pi)S l 


(4.3) 


we insert in the Green function the action, multiplied by a factor of e. Note that derivatives 
with respect to the couplings gi can be expressed as renormalized insertions of operators 
present in the action density 1 , 

-^7 exp (-S) = j d d zN{Ci(z)\ = J d d z^ (Sij + A ( z ) > ( 4 - 4 ) 

and then let us consider also the effect of the derivative 

d ^ d 
9 dg - d ( .h 

we obtain 


^<0|JV[Gi]^i...^ n |0) 



hence by taking the sum with the appropriate ~ e factor for the second term it is possible 
to cancel the “classical” contribution, that is the change in the Green function due solely to 
the measure. One has therefore 


d n d \ 
d/i dg) 

d 


= 0 


( 0 \N [Gi](j)i ... (f> n \ 0 ) = 
0 


2 eg— (EiMt j) Oj<j >i 



NG.j: 2e[ Pj - g ^ 


Sj<t>\ 



(4.6) 


1 in passim, observe that the action density is not a renormalized operator! 
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For simplicity, let us limit to minimal subtraction 2 , and have then 
L 

£P~ 


L p= 1 


Eq. 4.6 becomes 


(/^ + 2 eg-^j <0 | N [G^ ... cf> n \ 0} = 

J7 1 

= 2eg——~(0 \N [Gi\(t>\... 4> n \ 0) 



2eE i Y j g^ L O j (t>i...* 
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E 2 E 

j L,p= 2 
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£ p-\ S](p\ • • • 4>n 



where the “pre” /3 functions 


(4.7) 


(4.8) 


are introduced, the usual definitions of /3 and 7 functions will be recovered only when we 
shall specify the interaction part of the action density. 

By comparing with Eq. 4.4 we can modify the term multiplying the (3 in order to interpret 
it as a derivative with respect to a coupling: 


3 ft ^ln /?. 

g— + (pi + 2egp —— 2 eg ———(0 N [Gi]cf> 1 • • • 4> n \ 0) = 


dg 

= (0 
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d Mi i 
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00 , L 


2 E 1- 


L,p=2 


d \ P, 


L,p 


dg J e p 


-1 


-A— 

gi 


(4.9) 


We can now impose that the chosen subtractions make the Green functions finite: as the 
derivatives with respect to g and {g} do not introduce spurious divergences, we can say that 

2 Modified minimal schemes, like the MS or the G scheme, can be interpreted as a redefinition of the g 
scale and are therefore also minimal, in the sense that the finite subtractions needed to enforce them do not 
modify our formulas. 
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the l.h.s. in Eq. 4.9 is finite, and so must be for the r.h.s. First of all, in absence of the 
operator insertion, we should have then that 


r.h.s 







(4.10) 


is finite, and therefore each term multiplying the different 
follows the noteworthy relation 


pole parts must be zero. 


It 


2 ■LPj* 1 - 
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EaD 
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1 — > gi —— Loop counter operator . 

l d 9i 


(4.11) 


It will give useful relations between poles at different orders in L. The relation in Eq. 4.11 
allows to rewrite Eq. 4.9 as 


^ + (P l + 2e91 ) 


d_ 

dgi\ 


= \N [Gi](/)i 


9l 


(0 \N [Gj\(f>i 

<t>n\ 0 ) 


• 4>n\ 0) = 
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Ei (A + 2 e 9l ) ■. 
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■ 4*^ 


(4.12) 


Now we can impose that also the Green function with the insertion of the renormalized 
operator Gi is finite in the ee>0 limit. Taking into account that, by definition of minimal 
subtraction we have 

+ . ( 4 - 13 ) 

L,p 


we can rewrite with a few algebra Eq. 4.9 as 3 


+ ft + 2cm ) 

In /? 

= A—5—^(0 \N [Gi\4>i ...4>n\ 0 ) 


(0 | AT [Gi\4>i... <t> n \ 0) = 


dgi 


- 0 


E i '522LMZ’/N[G i ]<f> 1 ...<f> r , 


3 Observe that in presence of the operator insertion, which modifies the number of vertices, the loop counting 
operator is L = -g-§^ 
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(4.14) 


by defining 


7*,j 


= V 2 LM/y — - A 


L,lP _ 7 d In Ej ^ 
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we can rewrite the right hand side of Eq. 4.14 as 


(4.15) 
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and the last term has to be zero Vp, hence the recursive relation 
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In summary we have 


(4.16) 


(4.17) 
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(4.18) 


this form of the RG equation is written directly in terms of the “poles” of the Feynman 
graphs (the P and M matrices). 

It will be now useful to connect to the common notions of (5 and 7 functions. 


4.1.2 The conventional form 

Recall that in Eq. 4.1 we have introduced a coupling for each term in the action density, 
including for instance the kinetic term. This means that if we now want to set to unity 
the normalization of this term, we need to reinterpret the derivative with respect to this 
“eliminated” coupling: the usual way is to define a field anomalous dimension. 
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Following the constructive approach let us introduce the field anomalous dimension and 
proceed to its definition: for simplicity we limit the discussion to a single field </>, as the 
generalization is trivial. We rewrite Eq. 4.18 as 
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' , 7i + 


{pi + 2 £gi{ 


d n 
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7</> 77 j 


<0 \N [Gi\4>\ . ..(f> n \ 0) = 
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(4.19) 


and note that the n factor can be obtained by using the field counting operator 4 

5S 


h= d xN 


4>{x) 


$<t>(x) \ 

apart a correction due to the presence of the composite operator: one has 
n(0 | IV [0\(j)i... <f) n \ 0) = (0 | N [( 0}h<t>i... (/> n \ 0) - n(o 0-^lV [0\(j)i ...(j> n 
hence the r.h.s. of Eq. 4.19 becomes 
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(4.20) 


where we label with v l the number of fields in the composite operator Gi, and with pi the 
same quantity for a term O of the lagrangian density. Let it be g^ the coupling associated 
with the kinetic term d^cpd^cj), and define 7 ^ such that 

h = - 7 <t>g<t> {y 4, = 2 ) ; 

we obtain from Eq. 4.19, after the rearrangements needed to reassemble the normal prod¬ 
ucts, 

4 Note that the derivation is particularly simple in dimensional regularization, thanks to the simple form of 
the “action principle” 
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Figure 4.1: A contribution to Mj t ] 


d d n 

[ [ % + {Pl + 2£9l) Wi + 


= ^Ed 0 


d 


Pigi ~dgi + 1,1 ~ V3 • • • & 


<0 | N [Gi] </>!... </> n I 0) = 

o' 


7* 


<91n Ej 


7 « = 7 ‘ 7 + 1 R - a gi 
d 

Pi = Pi + i<t>—gi ■ 


s.. 


h3 > 


(4.21) 


The right hand side cannot be written in renormalized form and in fact is identically zero: 
indeed, let us consider a Feynman diagram with the insertion of the operator O* contributing 
to Mij: it has v J external lines, as in Fig. 4.1. 

Applying the operator p l gi-JE, it counts with a positive sign the number of lines connected 
to each vertex of the blob, and with a negative sign and a factor of 2 the internal propagators 

p ‘ s, wr ~ 21 + E p " • 

But one has also 
d + £ p v = 21 + d 


hence the identity 
P l 9l-^ + v l - d = 0 


(4.22) 


results. 
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4.1.3 Reference Formulas 

We have finally written the RG equations for 1PI Green functions (simply changing the sign 
of the term referring to the field anomalous dimensions ) as follows 




Si,j + 7 i,j \ (0 | N [Gj](j) 1 ...(f>n\ 0) 1PI = 0 


(4.23) 


and we have found the form of the (3 and 7 functions 
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(4.24) 


moreover, a set of recursive relations for pole parts have been found, and using the topo¬ 
logical relation in Eq. 4.22, they can be rewritten as 




£ 2iF L, r+ l 


E 2LM u M sf 
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g « - f 7) X ’; 


(4.25) 


these relations 5 can be expanded in h and give therefore useful checks on the computations: 
every pole in a diagram at order L in Ti. apart the simple pole, can be predicted in terms of 
the poles of graphs with ll < L and of the /3 and 7 functions. This is exactly equivalent to 
say, as we shall see, that order by order in perturbation theory the only new information, in 
the Renormalization Group Improved perturbation theory, comes from the simple log’s. 


4.1.4 Finite renormalizations 


Let us now consider the possibility of finite renormalizations, that is, we do not simply 
subtract poles but also finite parts. This can be useful in order to impose renormalization 
conditions, for instance on the matrix elements of operators. For definiteness, we shall suppose 
to have already regularized some operator basis {N [G]}, and we can therefore impose further 
renormalization conditions with a mixing matrix of finite parts. 

5 In proving the second one it can be useful to observe that in a manner similar to the proof of Eq. 4.22 one 
can show that 
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We introduce a matrix F and define the renormalized operators as 

R[G i ]=F iJ N[G j ]; (4.26) 

the matrix F is assumed to give finite corrections to counterterms, hence it possesses a 
loop expansion of the form: 

OO 

FiJ = kj + E F Zi • (4-27) 

L =1 

It is simple to find the modification to the RG equation Eq. 4.24 induced by finite renor¬ 
malizations: the RG equation can be written schematically as 
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and therefore the equation for the renormalized operators is 


FVF~ l + F^F~ l 


N 


R, 


= 0 . 


(4.28) 


(4.29) 


It is worth noting that the only source of + dependence is in the action, Eq. 4.1, hence 
the /i^ operator commutes with F. The derivative with respect to the couplings, instead, 
acts on F. and as we assume that this matrix is defined in order to give finite renormalization 
to counterterms, it has the same properties under this derivation as the matrix M: in other 
words the coupling content is the same, and it follows that (see note 3) 

g ®F = -LF. 

ogi 

By using also the identity 
FVF~ l = - (VF) F- 1 


the modified RG equation results 


V5 + f 
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Gi 


= 0 


7 = 


F^fF- 1 - (vF^j 
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-l 


dF 


= F* i F-L-\t3 l —-2eLF\F 


p-i 


(4.30) 


the effect of the finite renormalization is a modification of 7 in order to account for the 
evolution of the finite renormalization matrix, plus a “rotation” of the original 7 matrix. 
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4.1.5 Treatment of evanescent operators 

It is well known how in dimensional regularization the arbitrariness in the d dimensional ex¬ 
tension of 4 dimensional mathematical objects reflects itself in the appearance of the so-called 
“evanescent” operators. While we refer to the literature, and to the examples coming from 
practical computations (Chap. 5), for details on their treatment, in particular on the reduction 
formulas which are the equivalent of the Zimmermann formulas in the BPHZ formulation, we 
want here to discuss how evanescent operators modify the d—> 4 limit of the RG equation. 

Suppose for definiteness to be interested in the evolution of the matrix element of some 

operator N and that in the renormalization process it mixes with a set {N [R]} of relevant 

operators having a non zero limit in 4 dimensions, and a set {N [E] } of evanescent operators. 
The RG equation Eq. 4.23 are coupled, so in general one shall have 

VN R + rN R + rN E = 0 

VN E +7 e,eN E = 0 ; (4-31) 

now the set of evanescent operators is in many cases unbounded, that is, order by order 
in perturbation theory new operators appear; this is not the ruin of perturbation theory for 
in the d —» 4 limit the matrix elements of evanescent operators can be rewritten as local 
contributions, which can be therefore accounted for as finite renormalizations of the relevant 
operators. The well known result is that a reduction matrix exists, defined in perturbation 
theory, such that 

N[E] =r E , R N[R] ; (4.32) 

the matrix f is zero at the classical level. 

It is therefore immediate to decouple the evanescent operators in Eq. 4.31, by applying 
the reduction formula to the first relation in Eq. 4.31: 

vs + 7 '] n\r] =0 

l' = 7r, R T 7r, e^e,r (4.33) 

the second relation instead gives a consistency check: 

[(Pr E , r) - te, r (7R, R + 7R, e^e, r) + 7 e, e?e, r]R = 0 . (4.34) 

It is worth noting that the same results can be obtained with a slightly different procedure, 
namely by defining a non minimal scheme in which the evanescent operators do not contribute 
to the matrix elements, as their contribution is canceled by finite renormalizations. In practice 
this amounts to redefine the normal product 

N \e\ -> N \e\ - t e ,rN [r] (4.35) 

in such a manner to subtract not only poles proportional to other evanescent operators, 
but also finite parts proportional to relevant ones. 
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Let us find the resulting modification in the anomalous dimension matrix, as in sec¬ 
tion 4.1.4, by defining 

F = 5 - te,r , 

and then using Eq. 4.30 it results 

•j 1 = — r') 7 + rj + (Dr) (d + f^j ; (4.36) 

writing Eq. 4.36 in “components”, one easily finds that 

7R, R = 7R, R + 7R, ETE, R 

7 R, E = 7 R,E 

7 E, R = “PE, R (7 R, R + 7R, EPE, r) + 7E, EPS, R 

7e,e = 7e,e - 7’e,r7r,e • (4.37) 

A few comments are in order 

• the first relation shows that the “relevant” part of the anomalous dimension matrix is 
modified exactly as in Eq. 4.33; 

• the anomalous dimension mixing of relevant with evanescent is immaterial, as matrix 
elements of evanescent operators have been set to zero, and is unaltered; 

• we have written without comments in Eq. 4.31 that there is no mixing of evanescent 
operators with relevant ones, and consequently 7 e,r = 0. This is true in the minimal 
subtraction scheme, and if one takes into account the consistency condition in Eq. 4.34, 
one easily proves that this remains true, r = ^ otherwise, the matrix elements of 
evanescent operators would get non zero contributions from the relevant ones, vanifying 
our renormalization prescription, Eq. 4.35; 

• the last relation shows a modification of the flow of evanescent operators, again imma¬ 
terial in the d —» 4 limit. 

In summary, we have shown explicitly how two different ways to deal with evanescent 
operators, namely reducing them at the level of the RG equation or in the course of the 
renormalization process leads to the same results for the anomalous dimension matrix. 

We shall give explicit examples in Chap. 5. 

4.2 The RG improved perturbation theory 

In Chapter 3 we have seen how the Effective Hamiltonian resulting from the OPE defines the 
amplitude for a process i —> / as (we shall drop the 1PI suffix) 

= ( / K) complete 

« (■£? w)I w • 


(4.38) 




4.2 The RG improved perturbation theory 


38 


We have been able to compute the coefficients Cj for ji ~ M\y- Now we want to change 
the n scale in order to have /i~ m;, so to avoid large logarithms in the matrix element. 

We known that the fi dependence on the l.h.s. is equal to the one on the r.h.s. hence we 
can write down the identity 

( V + 7i + 7/) (/K)complete = 0 • ( 4 ‘ 39 ) 

On the r.h.s, the same identity holds, dropping unessential constants 

+ -* + V) Q (jfc, £;) (/ |JV [0 3 ]| i)*. M (4.40) 

but the derivative V acts both on the coefficients C and on the matrix element: we already 
know that a change in fi in the matrix element is compensated by a finite renormalization of 
the couplings and of the operators, that is, it holds the identity 

{V + 7* + 7/) i + 7] (/ |lV [o] | *) (4.41) 

hence the scale dependence of the coefficients results 

VC = 7 T C . (4.42) 

It is worth noting that there is no dependence on the external states, otherwise the effective 
hamiltonian framework would be meaningless. The Eqs. 4.41 and 4.42 are a mere tautology: 
any redefinition of the “normal product” N O can be compensated by a redefinition of the 

coefficient vector C, so to leave the physical amplitude invariant. 
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where one easily recognizes that the /? function is the Jacobian of the dimensional trans¬ 
mutation, and the P means “path-ordering”. 

In the LL approximation, a considerable simplification is possible. The coupling a factor¬ 
izes and the AD matrix can be diagonalized with a /j independent orthogonal transformation 
R 
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hence the path ordering is no longer needed and one immediately obtains 
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The evolution of the coefficients results 

We can observe that at leading order 
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(4.50) 


all the terms of the form a In have been summed, as expected. 

At the next to leading order, logarithms of the form a n In” -1 are summed, and so on. 
Detailed discussions of the NLL approximation can be found in a number of works (see for 
instance [22, 46, 77], and references therein). It is worth reminding that this successive ap¬ 
proximations are not directly connected with the loop expansion. For instance we shall see 
that the two-loop determination of the ADM in the b —» sj process allows only the resum¬ 
mation of the leading logarithms. This is due to a cancellation, at one loop one has no 
contribution to the amplitude of the form in Eq. 4.50. 


4.2.2 Scheme and scale dependence 


Physics is determined by the amplitude and therefore schematically by the combination 


C ■ N 


O 


therefore any scheme dependence implied by a different regularization and renor¬ 
malization scheme for the composite operators must be compensated by an appropriate re¬ 
definition of the coefficients. 

The scale (//) dependence itself is the reflection of the existence of a family of subtraction 
schemes, at different renormalization points, connected by the RG equation, whose coefficients 
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we are able to determine only approximately. In other words a scale dependence is introduced 
when we relate two subtraction schemes belonging to the same family, as the one used to 
compute the coefficients of the effective theory at the large mass scale Af to the one appropriate 
to evaluate the matrix elements at the lower scale ru¬ 
in fact there are three sources of approximation: the initial conditions of the evolution, 
at the higher scale Af, the coefficients of the differential equation, that is the (3 function(s) 
and the anomalous dimension matrix, and finally the matrix elements at the lower scale m. 

This scale dependence has nothing to do with the non-perturbative effects 6 but it is 
intrinsic of the perturbative expansion and of the reorganization of the perturbative series 
implied by the RG resummation. Assume for simplicity to have only one operator N \0\ 
contributing the some process i —> /, and to be able to sum the series 

{a n In”} , jct n In" -1 j , ,..{a n ln n - p } , (4.51) 

having computed the ADM up to the p + 1 order in the a s expansion. 

Then a consistent determination of the amplitude 

A(i-*f) = C (m) (f \N [O] | i) (m) (4.52) 

at the lower scale requires the knowledge of the initial conditions of the differential 
Eq. 4.44, C (Af, a) up to p order, and analogously for the matrix element (f \ N [0]| i) ( m ). 
For instance, terms aP from the matrix element, when multiplied by the terms of the series 
{a n ln n } in the coefficient, generate contributions of the form a n+p ln n . 

On the other hand terms of the form a n ln n ~ p_1 are not summed, in this approximation: 
this left-out logarithms are the source of the residual scale dependence. This dependence can 
be sized varying the lower scale around the typical scale of the i —* f process, for instance in 
the b —> sj case around the scale nib, and can be reduced, but never eliminated, only through 
higher order computations. 

4.2.3 Sizing the scale dependence 

In order to estimate the scale dependence it is necessary to consider the amplitude at the next 
to leading accuracy. For simplicity we shall limit ourselves to the case of a single operator in 
QCD (see for instance [46], and [73] for the general case), therefore the amplitude is 

C (Af, p, a) (f |N [0]| i) (m, p, a) (4.53) 

where we have made explicit the dependence on the coupling a, the two scales Af, m 
and the renormalization point p. If all the intermediate computations were exact, the p 
dependence of the matrix element would be canceled by the one of the coefficient. 

Let us assume that having p ~ m guarantees that no large logarithms appear in the 
expression for the matrix element: it means that it is possible to compute the p dependence of 
the matrix element accurately, even not using the resummation techniques. The p dependence 
is a purely short distance effect 7 , which means that a perturbative evaluation of the matrix 
element should suffice in checking the extent of cancellation of the p dependence. 

6 A further source of approximation, especially at the lower scale 

' Barring infrared complications, which in any case would be coped with including soft bremsstrahlung. 
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The strategy is as follows: we can write down the explicit expression for the Wilson 
coefficient, in dependence on M and fi, using the Renormalization Group, and check the 
cancellation of the // dependence against the matrix element, computed at the same accuracy 
in perturbation theory. This will give an idea of the implied accuracy. 

Given the NLL expression for the (3 function and anomalous dimension, 


(3qcd 
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(4.54) 


we easily write down the expression for the RG evolution operator 
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C (M, n, a) = U (n, M) C (M, M, a) 


(4.55) 


(4.56) 


and at the same accuracy we need to specify the Wilson coefficient leading a dependence: 

C (M, M, a) = (1 + • (4-57) 

Now the n dependence of the matrix element is easily obtained by expanding the RG 
equation to the appropriate order 


(f\N[0}\i)(fi) = U (m, n) (/ \N [0]| i) (m) 
U (m, n) = 
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i + — 7 °ln( 

VAt/. 


(4.58) 


Finally it is possible to write down the full expression for the amplitude, at NLL accuracy 8 
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8 The expression for the running coupling at NLL accuracy is 
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One easily recognizes that the dominant n dependence in the leading term is canceled by 
the logarithmic term coming from the evolution of the matrix element, as it should be. 

It is also important to note that in general it may be possible, given a NLL computation, 
to find an appropriate range in the scale fi which minimizes the size of the NLL corrections 
with respect to the LL ones: this is what one usally refers to as “setting the scale” of the 
leading computation by going to NL order. This is equivalent to say that for reasonable 
values of the NL coefficients in the 7 function, corresponding /x values can be found in order 
to reproduce the NL results with a LL expression. This somewhat “empirical” result helps 
to deduce the error implicit in the LL approximation by varying the /j scale in a reasonable 
range. The range itself however can be sized only explicitly checking the correspondence with 
reasonable values of 71 . 




Chapter 5 


QCD corrections to the b —> sj 
process 


We are now in position to discuss in full details the short distance QCD corrections to the 
b —> S 7 decay. This process has been studied by many authors and the LO corrections are by 
now well established. Nevertheless it is instructive to make a short review of the researches 
in this field. 

5.1 Review of the existing results 

The great interest of the b —» sy decay has been stimulated by the growing limits on the 
top mass. Indeed, the initial analyses based on the assumption of a large rrif had found a 
substantial enhancement of the decay rate due to QCD perturbative corrections; however 
different results had been obtained by groups working with different regularization schemes 
for 75 . 

The group of Grinstein et al. [37, 47], using the so called Naive Dimensional Regular¬ 
ization scheme (NDR), based on a fully anticommuting 75 even in d dimensions, found an 
enhancement of about a factor of 4, while the enhancement was limited to a factor of about 
2 in the work of Grigianis et al. [38], based on the so called Dimensional REDuction scheme 
(DRED) [15]. 

We stress that both the NDR and DRED schemes are known to be internally unconsistent, 
that is, to give wrong results in definite computations. For instance the NDR scheme does 
not reproduce the Adler anomaly, while the DRED scheme is known to fail at three loops (in 
4 dimensions) [18, 19]. 

It is quite easy to show the internal unconsistency of the NDR scheme: consider the 
following trace 

T^upa = TV Yl^al^lulplala] (5-1) 

and assume to first contract the a indices without exploiting the cyclicity of the trace and 
the assumed anticommuting property of 75 with the other 7 matrices: the result is 

Tpupa = {d- 8) Tr [757^7p7<t] (5.2) 
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on the other hand, bringing the 7 Q matrices close, one obtains 
T^pa = ~dTr [ 757 //)V 7 p 7 cr] • (5.3) 

Clearly the two results show that it is not possible to have a fully anticommuting 75 and 
consistently define its trace with 4 7 matrices. A possible way out is to maintain 75 a 4 
dimensional object, as in the t’Hooft-Veltman scheme (HV). 

This result does not mean that NDR (and DRED) schemes cannot be applied in many 
practical cases: as long as there is no need to give an explicit definition of the trace in 
Eq. 5.3, the NDR regularization will give unambiguous results. A caveat, all the subtleties of 
evanescent operators will still be present, see for instance [57, 76]. In fact various authors [22, 
73, 77] have been able to apply NDR, DRED, HV schemes successfully to the QCD corrections 
to weak hamiltonians, obtaining consistent results. 

As the NDR scheme does not appear to cause unconsistencies when applied to b —■> sj 
decay at the leading log approximation, we confirmed in [57] the results of Grinstein et al. 

All these works were performed in a reduced basis, not including the mixing with a certain 
class of 4 fermion operators. The first author to extend the computation to include all the 
mixings relevant in the LO approximation has been Misiak [58], using the NDR scheme. 
The resulting corrections to the previous NDR computations are small, showing that the 
approximation used in [37] was good. 

Nevertheless the use of a complete basis has led to understand the origin of the dis¬ 
agreement between the pioneering works. Ciuchini et al. [76] have shown that the operators 
completing the basis contribute to the b —> sy amplitude, through one-loop graphs, at zeroth 
order in QCD. This contribution is scheme dependent (for instance, in the HV scheme is 
absent), and compensates the scheme dependence of the AD matrix. 

We shall fully describe the re-evaluation of the LO QCD corrections [102, 103], based 
on the scheme recently used by Curci et al. [74], a method we shall call from now on the 
“symmetrized” scheme. In our opinion this technique greatly simplifies the calculation and 
will be useful in the NLL three loop computation. 

5.2 The symmetrized scheme 

In the Chapter 3 we have derived the Effective Hamiltonian for the rare B decays: here we 
add the operators required by the QCD corrections 
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(5.4) 


In principle working off-shell they also appear non gauge invariant operators, which do not 
contribute to physics as they are proportional to motion equations. This has been avoided 
by using a background gauge fixing term for the gluon field, thus ensuring that the gauge 
invariant basis in Eq. 5.18 closes under renormalization. 

The coefficients different from zero have already been given in Eq. 3.13, and it is to be 
noted that they are evaluated in the NDR scheme: in the LO approximation the scheme 
dependence of the coefficients at the n — Mw scale is irrelevant, at least for the 6 - 4-57 
process. 

Two comments about the basis in Eq. 5.4 are needed: first of all this off-shell basis can 
be reduced, using the equations of motion, and we shall elaborate on this in the following. 
Secondly, we have seen that in the d dimensional formulation of the RG equation they are 
present also evanescent operators. In what follows we assume that a subtraction scheme is 
chosen in order to decouple them from the coefficient evolution: this will be explicitly shown 
where appropriate. 

Let us consider the general structure of the Effective Hamiltonian, leaving aside the field 
content 


n eS = ^cf w iv[i? i ] + ^c'^ w iv[L i ] 

i i 

+ E C^ HV N [(L ® L)-] + E C^ HV N [{L ® R\] ; (5.5) 

i i 

N [L, i?] stand for operators bilinear in fermion fields, like the magnetic momentum oper¬ 
ator, and N [L <S) Li], N[L® Ri ] stand for current-current operators. 

The symbol R and L reminds the presence of the chiral projectors Pl = |(1 ± 75 ). This 
effective hamiltonian is obtained from the “complete” theory, say the SU ( 2 ) x U{ 1) model: 
if we make in the complete theory the substitutions 


75 


75 = -75 
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(5.6) 
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the resulting effective hamiltonian (in the same regularization and renormalization scheme) 
is identical to the one in Eq. 5.4, except for the same substitutions in Eq. 5.6. 

But the QCD corrections are independent on the convention used for defining 75 , hence 
the RG evolution equation is invariant under the trasformation (5.6) 


' d ; 

M-r- + 7 
a/i 


/ N[L®L] \ 
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d : 
L% + 7 . 
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(5.7) 


In other words the anomalous dimension matrix 7 is the same for the two set of operator 
in Eq. 5.7 and for the purpose of the calculation of the AD matrix any linear combinations of 
the two sets is equivalent. It is convenient in particular to consider the RG evolution of the 
symmetric combination, 


' d ; 

V-T + 7 
a/i 


( ± (N[L® L\ + N [R®R]) 
| (N [L®R]+N[R<g) L}) 
V ±(N[R + L}) 


= 0 , 


(5.8) 


because it is possible to redehne these operators in order to make the 75 matrix disappear. 
This is trivial for the symmetrized magnetic momentum operator, while for the 4-fermion 
operators a change of the d dimensional extension is needed. 

A complete (infinite) basis for the Clifford algebra in d dimensions is given by the com¬ 
pletely antisymmetric products of 7 matrices [20, 30, 56, 74], 


7^ ^ — 7/11,M2,-Mn — I 5Z ( 1) P 7 /u 7£12 • • • 7/in • 
n • pen„ 

In 4 dimensions the structures 7 ^ survive only for n < 4, and one can write 


(5.9) 


^ [(L (g) L) + (R (g) R)] = ( 7 ^ < 8 > 7 (1 ^) + ^j(7^ ® 7®) 

^ [(L <g> R) + (R (g) L)\ = <g) 7 (1 ^) - ^(7^ < 8 ) 7 ®) • (5.10) 

These equations can be taken as definitions for the tensor products in d dimensions; in 

practice we change the d-dimensional extension of the “symmetric” operators, with the aim 

to simplify the evaluation of the AD matrix. 

The key advantage over the HV scheme is that there is no split between 4 dimensional and 
d — 4 dimensional objects (the indices in Eq. 5.9 are d dimensional). For example, working in 
the well definite basis of ^y( n ) 0 7 ^)^ structures, for arbitrary integer n, relevant operators, 
having n < 4, and evanescent ones having n > 5 are treated in a unified manner. 

We stress that the definitions in Eq. 5.9 render this scheme different from the HV one. 

In fact, working in the HV scheme, the effective hamiltonian in Eq. 5.4 can be cast in the 
form 


Keff = E cf v {n)N[ot\n) + c? v ^)N[o\~\n) 


(5.11) 
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where O £ and oj ^ are the symmetrized and antisymmetrized operators. The two classes 
of operators are not mixed by QCD corrections, and so the two terms of 7i e ff are separately 
RG invariant 

l‘^T, C . HV ^) N Pi +) M = ^Y.C, HV MNlot%)=0- (5.12) 

The redefinition of the symmetric part amounts to add of evanescent operators to 
the 0^ + \ defining the “symmetric” and extended operators 

0^=0\ +) + E^ . (5.13) 

Using the operators one obtains an AD matrix different from the one obtained with 
the operators, which are defined as in the HV scheme. The mismatch must be readsorbed 
in finite renormalizations, which ensure that we can obtain the AD matrix in the HV scheme, 
and therefore also the evolution of the antisymmetric part. 

This is possible as a renormalized evanescent operator can be expanded on a complete 
basis of relevant ones [27] with finite reduction coefficients r t j , 


N 


E, 


(+) 


*) = nj(f 


N 


O 


(+) 


For instance imposing the condition 


£c££)IV 0, (+) (M) = £cfV)iV 


(+)/ 


o 


(+) 


(m) 


(5.14) 


(5.15) 


one has 


c\ + \n) = (8 ij -r ji )C» v {v)- (5-16) 

The RG evolution of coefficients is guided by the anomalous dimension matrix 

governing the evolution of the symmetrized basis, and using the Eq. 5.16 it is simple to find 
the evolution of Cf v (n). 

Alternatively the Eq. 5.16 can be interpreted as a redefinition of the normal product N \0\ 
as a non-minimal subtraction procedure: we shall use the notation N' [O] for this prescription. 

The advantage of the method is that the coefficients are O (fi); hence the most cum¬ 
bersome part of the computation, the determination of the AD matrix, can be done in the 
symmetrized scheme, while the matching needed to write down the amplitude requires a 
computation at one loop order less. 

In the following we shall use the following symmetrized and extended basis 


(+) = 

v/2 


7-/vw = 

n eS — 


v t * a v tb 


£3 


(+)q(+i 


(5.17) 
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where the operators are defined 

of = f^spppb 
of = <=§$* 57* bD^ 

of = flf m b s b 
of = f$s^T A b(D,G, w ) A 

of = ^(s a p n) C p) ® (cp 7 (n) M 
Olli = ^(Sa7 (n) fta) ® J2q(Ql3l in) qi3) 


Of = S {P, 0>} & 

Oi +) = -( 3 \?m b sppb 

Of = §$*{A TA G A v °»v) b 
Of = fjf-m, b sT A G A u a^b 
Of n = hff n)c 0) ® (C/37 (n) M 
Oljl = ® Egfer^a) 


(5.18) 


and the non-zero coefficients are normalized as follows, with respect to the ones in Eq. 3.13: 
of = 2 Ci 1 = 1,... 8 Cj, +\ = Cf,J } 3 = Cm • (5.19) 

5.2.1 On-shell and off-shell formulation 

To make contact with the work of other authors [76, 58], we shall also use an on-shell basis, 
obtained from Eq. 5.4 by exploiting the motion equations. The reason for computing QCD 
corrections off-shell is that the redundancy of the computation provides further checks, and 
the intermediate results can be saved in view of the planned NLO computation: nevertheless, 
given the AD matrix, it will be convenient to reduce the results to the following on-shell basis. 


Hi 

Qi 

Q-2 

Q 3 ,5 

Qa, 6 

Qi 

Q8 


t r-k t r \ A /-yon— shell r\ 

Qi 

( sc) v _ A ( cb) v _ A (1 <g> 1) F 

( SC)y_ A ( Cb)y_ A (1 <g> 1) 

( sb) v _A !) 

q=u,d,s,c,b 


{sb) v _ A ^ ( < m)v±a (1 ® 1 )e 

q=u,d,s,c,b 

(— icQd) _ __ jj , ^ 

, \2 [iv (V+A)FqLvb 

(47r) 

{-igs) — a , 

, . .9 rn b s<J fj,n (V+A)'^ r iu'b ■ 

(4tt) 


(5.20) 


It is worth recalling the general properties of an anomalous dimension matrix 7 for an 
operator basis not reduced by motion equations: following [5, 6 , 7], it has the block form 


relevant 

motion equations 


relevant 


X 

0 


motion equations 


Y 

Z 
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and we shall explicitly check this result in the computation. Please note that QCD cor¬ 
rections background gauge fixing term has been used for the gluon field, thus ensuring that 
even off-shell there is no need to introduce non gauge invariant operators, and that the basis 
in Eq. 5.18 closes under renormalization. 

The symmetrized and extended version of the basis in Eq. 5.20 will be given by the 
following combinations of the operators appearing in Eq. 5.18 



5.3 The computation 

5.3.1 One loop results 

We list in Fig. 5.1 and in Fig. 5.2 the general structure of the graphs needed to renormalize off- 
shell the operators 8 . Here and in the following figures the external zigzag line addresses 
a gluon or a photon. Note that the graphs in Fig. 5.2 are proportional to aj, but as pointed 
by Misiak [58] they are required even at leading order since they give rise to the mixing of 
the operator O 1 * * * * 6 7 8 with the operators n) ^ 12 n y 

The computation of the Feynman graphs in Fig. 5.1 results in the {1... 8} x {!-•- 8} 
sector of the one-loop anomalous dimension matrix qh), 

1 2 3 4 5 6 7 8 

1 /-2 c F z 0 0 6 C F 0 0 0 0 

2000 0 8 C F 0 0 0 

3000 0 00 0 0 

cn s 4 0 0 0 6 Cf-2 C f £ 0 0 0 0 

47t 50 00 0 8 Cp 0 0 0 

6 -12 C F 4 C F 0 0 4C F 4C F 0 4(3 C F -C A ) 

7 2 C F 0 0 0 0 | Cp+3 Ca 0 

8 \ 0 0 0 — 12Cf 8C f 0 0 4(4 C f —Ca) 
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Figure 5.2: Mixing of operators Of 8 with 0( t L1 _ 12 n \ 


$MW6~(T6'(r6Mm6l 
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Figure 5.3: Mixing among 4-fermion operators. 



Figure 5.4: Mixing of 4-fermion operators with 0+ 8 operators. 


while the graphs in Fig. 5.2 connect the operator Oj to four fermion operators, resulting 
in 



(H,l) 

— 5+5 C 2 a —10 C A Cp 
2 


(11,3) (12,1) 

9 C\ -30 C' A C F +24 C% n 

2 U 


(12,3) 

2iC F -9C A 

2 


(5.23) 


The renormalization of four fermion operators, off-shell, results from the graphs in Fig. 5.3 
and in Fig. 5.4. 

The self mixing of four fermion operators at one-loop, in Fig. 5.3, connects Dirac structures 
with An = +2, 0, —2 


Os_ (9, n) 
(10, n) 


(9, n — 2) 

(6—n) (n-5) 
C A 

(n— 6) (n— 5) 
2 


(9, n + 2) 

— ((1+77.) (2+n)) 

Ca 

(1+n) (2+n) 

2 


(9, n) 

2 C F ((1—n) (n-3)-0 
3 (—2+4n—n 2 ) 


(5.24) 
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(10,n-2) (10, n ) 

Oi s (9, n ) ( (n—6) (n—5) 0 

47T (10, n) y ( 4 C F - CA )( n - 6 ) ( n - 5 ) ( 3 C A - 4 C F ) (2-4n+n 2 )-4 C F (1+0 


(10, n + 2) 

(l+n) (2+n) 

(4 Cp — C A ) (l+n) (2+n) 
1 


The same mixings occur in the n) ^ 12 sector. 

The penguin graphs in Fig. 5.4 give rise to the mixing of the 4-fermion operators with 
Oil' 7 at the a° order 


3 7 

(9,n) / -^(-l) n c n (0) ^(-1)^(0) 

(10, n) -C A ^=p-(-l) n Cn(0) 0 

(11, n) C A %n f 6 nil + ^^(-irc n (0) ^^(-lfc n (0) 

(12, n) \ l^/^n.,1 + C , A 2( 3 ~ ! 2) (~ 1 ) n Cn(0) §7+/<5n,l 

The meaning of the symbol c n is explained in App. B together with the other definitions 
and useful formulas. We use the symbols rif = u + d, fif = d — 2u, with u and d being the 
number of up and down quark species active. 

It is to be noted that at leading order only the n = {1, 3, 5} values are needed, because 
the effective hamiltonian starts with n = {1, 3} and at one loop only the n = 5 evanescent 
arises. We find convenient to give also the results restricted to this set of n values 



(9,1) (9,3) (9,5) 
(9,1) / 0 0 

a, f (9,3) 0 -g 

47t (10,1) 3 3 0 

(10,3) 3 3 10 


( 10 , 1 ) 

0 


(10,3) 

6 


(10,5) 

0 


"UX + 6Cf 


-^- + 6C i f 0 

*C F 20C F -£ 


(5.26) 


3 7 

(9, !) / -I | 

(9) 3) ~ 3 3 

( 10 , 1 ) 0 

(10, 3) 0 

(11,1) |(1 + C A n f ) | 

(11, 3) 3 3 

(12, 1) § (C A + n f ) | n f 

(12, 3) V | C A 0 

It is well known that renormalized operators proportional to motion equations mix only 
between themselves and that their anomalous dimension matrix is not gauge independent. 
In fact, we can observe that operator N , proportional to the s motion equation, mixes 
only with itself and operator N [OJ]: this one mixes with itself only, and both operators have 
anomalous dimension matrix depending on £ [5]. The difference N \ Og" — Og" 1 is proportional 
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to a combination of the s and b motion equations and does not evolve at one-loop. Analo¬ 
gously the operator N Og" — O g~ mixes with itself and with N , N — O 5 " , N 0~l . 

Finally the combination 

^ + £(2 

is proportional to the equation of motion of the gluon. It is worth noting that the elim¬ 
ination of the operator N Oj in favor of the four fermion operators introduces a factor 
which combines with the o? s in Eq. 5.23 to give a result relevant for the leading order 
computation [58]. 

5.3.2 Two loop results 

The two-loop mixing of the four fermion operators with operators 8 can be obtained from 
the computation of the Feynman diagrams in Figs. 5.5, 5.6 and in Figs. 5.7, 5.8: the “closed” 
loop graphs are relevant only for the renormalization of operators 0 ^ n _ 12 ) n y 

The massive b quark is represented by the heavy-faced lines. The b propagator is expanded 
in series of mb up to the first order and the resulting massless graphs are regularized in the 
infrared region by the flow of the external momenta: the diagrams have to be evaluated with 
zero and one mass insertions. The loop integrals have been computed with the help of the 
algorithms developed by Chetyrkin et al. [21], implemented in the Mathematica [62] symbolic 
manipulation language. As in the case of the one loop computation, the results can be given 
for arbitrary n and are listed in Tabb. 5.1, 5.2, 5.3. We refer to App. B for symbols and 
definitions; let us just note that the symbol f n j results from traces involving elements of the 
{ 7 O 1 )} basis. 

We stress that having the results for arbitrary n will be useful for the NLO computation. 


5.3.3 Reduction of evanescent operators 

Before giving the results for n = 1,3, that are needed in this leading logarithmic computa¬ 
tion, we perform the reduction of the evanescent operators with n = 5. We follow the scheme 
detailed in Sec. 4.1.5, computating the graphs in Fig. 5.9, which allow to express the inser¬ 
tion of the evanescent operators in the Green functions in terms of the insertion of relevant 
operators, with coefficients addressed by the matrix f 


(9.5) 

(10.5) 
■ (11,5) 

(12.5) 



(5.28) 


For the present computation only columns 5, 8 are relevant, while the columns 3, 7 con¬ 
tribute through equations of motion to four fermion operators, 


(H,l) 

(9> 5) / 10 C A 

= a* (10,5) 0 

- 47T (11, 5) 5 ^ 

(12,5) V 0 


( 12 , 1 ) 



(5.29) 
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Figure 5.6: “Open” graphs for the renormalization of 12 n y Non abelian couplings. 


and are relevant only for the NLO computation. 

By combining the reduction coefficients in Eq. 5.28, the one loop anomalous dimension 
matrix in Eq. 5.26 and the results of the two loop computation listed in Tabb. 5.1, 5.2, 5.3, 
we are able to give the anomalous dimension matrix relevant at leading order, where a factor 
ft is understood. 

4 7T 



1 

2 

(9, 1) 

/ 2 C F 

3 

—44 C p 
9 

(9, 3) 

2 C F 

9 

—44 C p 
9 

(10, 1) 

0 

0 

(10, 3) 

0 

0 

(11, 1) 

4 Cf 

3 

20 Cf 

9 

(11, 3) 

4C F 

9 

20 Cf 

9 

(12, 1) 

Q 

to 

Si 

1 

CO 

S 

16C f (2n r : 


27 

(12, 3) 

l o 

4 Cf n f 
3 


—88 Cp 

8 Ca-^t^- 

—40 Ca Cf 
3 

4 C A (3-14 C F ) 

3 

-176 , 8 C A C F n f 

27 3 

944 Cf 
27 

40 C A C F | 8C F (95n f -102n f ) 
3 81 


32 Ca Cf 


+ 


l Cf rif 


4 

-2 C F \ 

—2 Cp 
0 
0 

-16 C F 
-24 C F 
4Cp (2nf—3nf'j 
3 

0 


5 6 


(9, 1) 

/ 

—8 Cp 


Ca i 10 C F 

2 ' 9 


\ 

(9, 3) 


-8 C F 


C A , 10 C F 

2 ' 9 



(10, 1) 


0 


1 



(10, 3) 


0 


1 



(11, 1) 


-4 C F 


C A + ^ 



(11, 3) 


-32 , 32 C A , 28 Cf 

3 Ca 3 1 3 

C A + 

*f* + 2(n f - 


1 

(12, 1) 


4 C a Cf 

0 , (16C F -18C A )(2n / - 
w + 27 

■ 3n /j 


(12, 3) 

V 

52 C A C F I 8 C F n f 

3 1 3 

2 + 

(4Cf—Ca) (3nj--2 
3 

n f) 

/ 
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Figure 5.8: “Closed” graphs for the renormalization of 0^ 1 _ 12 n y Non abelian couplings. 





Figure 5.9: Diagrams relevant at LO for the reduction of evanescent operators. 
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Table 5.2: Two loop AD matrix entries for n y resulting from the mixing with Of g 
7 entry 

( 11 , n), 1 2 Cf 3 ( n T 1} [(1 - n) c n ( 0 ) + (2 - n) c n ( 1 )] 

( 11 , n), 2 [/„, 3 ( 0 ) + | (2 - n) (-l) n c n ( 0 )' 

( 11 ,n),3 ^[8C A n f 6 nA + 4f n:3 (0) 

+| (-l) n (n-1) n (—1908 + 1124n — 213n 2 + 13 n 3 ) c n _i(0) 

+£ (-l) n (—58 + 230 n — 177n 2 + 33n 3 ) c n (0) 

+| (-l) n (n - 2) (-11 + 24n - 6 n 2 ) c n (l) 

(11, n), 4 [i(39- 68 n + 20n 2 + 3 (2-n)(-l) n ) c n (0) + 2 (n - 3) (n-1) c n (l) 

(11, n), 5 ^ 2 (n - 6 ) (n - 5 ) 2 (n - 3) (n - 1) nc„_ 2 (l) 

-± (45 + lOOn - 170n 2 + 56n 3 - 4n 4 

+2 (n — 2) (2n 2 -8n + 3) (-1)”) c n (0) 

+2 (n — 3) (n — 1) (—5 + 3n — n 2 ) c„( 1)] 

( 11 , n), 6 ^ 12n/5 n , 3 + | (4Cf - CU) / n , 3 ( 0 ) 

+( 9 C A ~ 8 C'f) (n - 2 ) c n ( 0 ) 

( 11 , n), 7 ^ [ 88 nf5 n $ + (4 C F - Ca) /n, 3 ( 0 ) 

+ ±(-l) n ( 2 C F — Ca) (n- 1 ) n (-1908 + 1124n - 213n 2 + 13n 3 ) c n _ 2 ( 0 ) 
+Bi(-l) n (3 C A (n — 4) (-3 + 46n - 24n 2 ) + 

+ 4 C F (-58 + 230n - 177n 2 + 33 n 3 )) c n (0) 

+±(-l) n (2-n) (3 C A (3 + 8 n - 2n 2 ) + 4<7 F (11 - 24n + 6 n 2 )) c n (l) 
( 11 , n), 8 ^ \{-Ca + 2 Cf) (n - 6 ) (n - 5) 2 (n - 3) (n - 1 ) nc„_ 2 (l) 

++ (CA (99 + 24 n — 125 n 2 + 40 n 3 — 2 n 4 ) + 

+ 2C f (—45 — 100n + 170n 2 — 56n 3 + 4n 4 )) c n (0) 

+|(n - 3) (n - 1) (Cau + 2CV (-5 + 3n - n 2 )) c n (l) + 24 rif <5 n , 3 

(—1)™ (n — 2 ) (8(7^ (—3 + 8 n — 2 n 2 ) + ( 7,4 (l — 16n + 4n 2 )) c n (0)j 


al. [76, 77] in the HV case, 


Qt 

= Qf—Ef = 

Qf 

= Qf Ef = 

Q 3/5 

= Qf/ 5 “ E 3/5 

Q4/6 

= O 4/6 _ ^ 4/6 

Qf 

= Qt = of 

Qt 

= Qt = of 

Qt 

= Qf-E+ = ~ 

Q10 

= Qto — Ef 0 = 


>( 9 , 1 ) + 0 ( 9 , 3 )-^ 

> 70 , 1 ) + 0 ( 10 , 3 ) ~ ^2 
= 0(1!, l) + O^ 3 ( - 77y 5 

= 0(i 2 , i) + 0 ( + 12; 3 ( - 


79 , 1 ) 

9 + 

- 710 , 1 ) 


i+ _ p+ 

710 , 3 ) ^10 • 


(5.30) 


The last two operators have been introduced to have an invertible relation between the 
bases Of and Qf ; as expected they decouple from the others in the RG evolution. The 
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Table 5.3: Two loop AD matrix entries for 0^ 2 n y resulting from the mixing with 0 2 j 
7 entry 

(12, n), 1 C F nf5 nt i 

(12, n), 2 — y CFnf5 n ^ + 4:CFn,f5 n) 3 

(12, n), 3 — gj Cf (102n/ — 81h/) 5 n , : i + ^ Cpnf5 ni3 

+ CAC 9 r n ( !~ 1)n t( n - !) n (-1908 + 1124 n - 213n 2 + 13n 3 ) c n _ 2 (0) 

+ (-108 + 148n — 83n 2 + 13n 3 ) c n (0)] 

(12, n), 4 -4C , Fn/<) n ,i 

(12, n), 5 |(n - 6) (ra - 5) 2 (n - 3) (ra - 1) ra c n _ 2 ( 1) 

+| (45 - 80n + 56n 2 - 12 n 3 ) c n (0) 

+ (?r — 3) {n — 1) (—2 — 5 ra + n 2 ) c n (l)] + 8 Cf fif 5(n , 3) 

(12, n), 6 / n ,3(0 ) + | (9 C*a — 8 Cf) ray <5 ni i + (4 Cf — C^) ray <f nj3 

(12, n), 7 §§ (9C^ — 2 Cf) nf5 n , 1 + ^ (4 Cf — 6(4) nf5 ni3 + £jf n , 3(0) 

ppp [l (n- 1) n (—1908 + 1124n — 213n 2 + 13n 3 ) c n _ 2 (0) 

+l(24 + 32n-35n 2 + 9n 3 ) c n (0) 

+2 (2 — n) (3 — 8 n + 2?r 2 ) c n (l) 

(12, n), 8 5 3 C ' 4 n/(5 n ,i + 2 (4Cf - Ca ) nf5 n ^ 

+ 37TI 5 ( n “ 6 ) ( n - 5 ) 2 ( n -3) (n - 1) ra. c n _ 2 (l) 

+1 (73 ra 2 - 18 - 62 n - 22 n 3 + 2 n 4 (-l) n (18 - 57ra + 36 ra 2 - 6 ra 3 )) c n (0) 
+5 (ra - 3) (n - 1) (-4 + lira - 3ra 2 ) c n (l)j 


presence of the evanescent operators takes into account the difference of the d-dimensional 
extensions. In the basis Qf we give the final results for the AD matrix in the following block 
form: 


a, / 7// 7/m \ 
47T y 6 7 mm J 


(5.31) 


The matrices 7 //, 7 mm , result from the one loop computation and are scheme independent; 
they guide the RG evolution respectively in the four fermion and in the magnetic momentum 
sectors, 


7/7 



7 f 8 C F 

8 y 8 Cf 
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6 
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0 
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2 
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2 

2 

C A 

3Ca 

3 

3Ca 

3 

0 

22 

22 

4 

4 

3Ca 

3 

3Ca 

3 

0 

2(9 C A -n f ) 

2(9 -n f C A ) 

2n f 

2n f 

3 C A 

3 C A 

ZCa 

3 

0 

0 

0 

6 

C A 

-6 

0 

2 n f 

2 n f 

2n f 

2(n / -18C F ) 

3C A 

3 

3 Ca 

3 


8 

0 ) 

16 Cf — 4Ca J 


5 


7mm 


(5.32) 
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The matrix 7 f m connecting the four fermion and the magnetic momentum sectors results 
from the two-loop computation and it is scheme dependent, even at leading order, as pointed 
out by Ciuchini et al. [76, 77] 


7 fm — 


V- 


0 

232 Cp 
9 

280 cy 

MCkCk + 4C £ 9 ( 27 % _ 4n/ ) 

104 Cp 
3 

40 Ca Cp 4 Cp 
3 9 


(Arif + 27 rif ) 


6 \ 

8 C A I 92 C F 
3 ' 9 

^ + 6n f 

-4 -*§*n f + Z*f*n f 

32Ca _ 104Cp_ 6n/ 


(5.33) 


5.3.5 Matrix elements and scheme independence 

The scheme dependence of the matrix 7 j m is compensated in the physical amplitude by the 
matrix elements of the four fermion operators. 

The “even” contribution to the amplitude for the b —► S 7 decay has the form 


< s'rlH+flb >= C 7 S (s 7 | N [Q+] I 6) o + jh C * (s 7 \N [Q+] | 6), 


n =3 


(5.34) 


where the subscripts 0 , 1 address the order in h of the matrix element, and an analogous 
formula holds for the b —> sg amplitude 1 . 

The one loop matrix elements of the operators 4 5 6 result from Feynman diagrams 
analogous to the ones listed in Fig. 5.9 and with a mass insertion. If we define our scheme 
preserving the naive four dimensional Fierz simmetry they give a zero result. In fact in this 
case the “open” diagram in Fig. 5.9 with a mass insertion is proportional to the “closed” one, 
which is set to zero by the trace. The symmetry can be preserved by using the HV scheme 
for 75 and maintaining the four dimensional definition for the current-current products [77]. 

In other regularization schemes this is no longer true and in general one obtains a local, 
finite contribution. This difference can be reinterpreted as a finite correction 2 to the coefficient 
Cf or alternatively to the operators. In our scheme the following relations hold 


s 7 


sg 


0 _ n, n 
®12 ,n 


b) = 


1 

(~36+64n—19n 2 ) c „(q) (1-n) (-3+n) c n (l) 

n! ' 6 n! 


36 


o + 

°T1, n 

f)+ 
°T2 ,n 



1 

Ca 



/ (—36+64n—19n 2 ) Crl (o) 
l 36 n! 


1 (1-n) (-3+n) Cn(l) 
6 n! 



(5.35) 


x The bar over the coefficients reminds us that they are obtained after applying motion equations: for 
instance C? = C| + C a 5 , = C a 6 +C a s . 

2 see, f.i. in [81], the so called “effective coefficients” formalism 
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Using the definitions in Eq. 5.30 we can write 



Z7i(s'y\Qj 

Zsi (sg I Os 


i = {1...6} 
i = { 1 ... 6 } 


where following [76] we define the vectors 


(5.36) 


Z7 

Z 8 


( 8 8Ca 2 2 Ca 

V ^ 3 5 3 3 ’ 3 




(5.37) 


We can perform the finite renormalizations needed to compare the results with the ones 
obtained in the HV scheme. We match the renormalization schemes by defining a non-minimal 
subtraction 


N' 


Ql 


5F l3 


{ l + 5P ) l3 N 




-Zji j = {7, 8 } 
0 j = {!,... 6 } 


(5.38) 


which sets to zero the matrix elements in Eq. 5.36, as they are in the HV scheme. 
According to the formula in Eq. 4.30 the AD matrix is modified as follows 

f = (l + SF^(l + AF )” 1 . (5.39) 

Using the results in Eq. 5.37, one easily obtains the AD matrix in the HV scheme 


7m/ 


1 

2 

3 

4 

5 

6 


7 

( 0 

208 C F 

9 

232 C F 
9 

^n f + 12C F fi f 

— 16 Cp 

V n f ~ 12 Cf ' n f 


8 

6 \ 
— AC A 

-8C A + 6n f 

12 +(iif 7 _ 4Ca ) n f 

AC a — 16 Cf — 6 nf 
- 8+(2 C A -m^)n,J 


(5.40) 


which coincides with the result in [76]. 
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5.4 Effective Hamiltonian at /i ~ 


We have obtained the expression for the AD matrix in the HY scheme: 


a. 


7hv = 
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0 
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3 
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0 

0 

0 

0 

0 

0 

32 

3 

28 

3 
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(5.41) 


where the values rif = 5, hj = —1 have been inserted, as appropriate to the 5 quark 
theory below M\v- 

With the help of Mathematical 62] it is easy to find numerically eigenvectors and eigenvalues 
of the transpose of the matrix in Eq. 5.41 




(5.42) 
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-0.02363 

0.04644 

-0.02528 

0.1849 

-0.3489 
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0 

0 
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0.7071 
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0 

0 
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0 
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Then we can determine the evolution of the coefficients at leading logarithmic order, on 
the basis of Eq. 4.49 


C (mb) = R 1 


/ a s (m, b ) 
Vets (M w ) 


— D 

2 / 3 0 


■R-C(Mw) • 


(5.45) 
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The numerical result for the coefficient C ™ 1 shc11 of the magnetic momentum operator, is 
at n = mb 


r = 

C% n - sheU (m b ) = 


+ 


«s (m&) 

(%) 

r -0.69565 c on-shell + g ^-0.69565 _ r -0.6087^ ^on-shell 

^6.52786 ?’ -0 ' 69565 - 13.797 r “°- 608696 + 3.896 14 r -0 - 40862 

+ 2.57143 r “°- 26087 + 0.034332 r -0 14565 + 0.22793 r 0422989 
+ 0.428571 r °- 521739 + 0.111299 r 0 - 899395 ) C% n ~ shel1 (M w ) . (5.46) 


The values 3 of the on-shell coefficients can be obtained from Eq. 3.13, taking into account 
the change in the normalization of the operators 


Con-shell ( Mw ) = x 

= 2 (C 2 (M W ) + C 5 (M W )) = - 71 - 3 f/ 31 , 2) In, 

4 (x — 1) 2 (x — 1) 

C° n - shell (A%) = 2 (C 6 (A%) + C 8 (+%)) = ' 1 1 + lna; - 

4 (x — 1) 2 (x — 1) 

Following Buras et al. [81] we shall normalize the QCD coupling to the value determined 
at LEP, 


(M z = 91.175GeV) ~ 0.12 


(5.47) 


using consistently the leading order approximation for the running coupling, 

/ x Oi s {Mz) 

as W ~ 1 a,(Mz) 11 C A —2rif , f ' 

1 "T 2t r 3 111 V M z ) 

With this choice the ratio r has, still for /j, = mb, the value 


(5.48) 


r ~ 1.7 (5.49) 

with the choices mb = 5GeV, Mw = 80.14, and the formula for the coefficient ( 7° n_she11 
becomes 

c on- she Il ( mb) _ Q 97 c on-shcll + Q gg ^on-shell _ Q.26 ^on-shell . (5. 50 ) 

as anticipated the multiplicative renormalization tends to suppress the ratio. 

Considering the range of variation of the coefficients for in the 130 — 180GeV interval, 

c on-sh e ii {Mw ) = o.95 - 1.20 

C™- sheU {M w ) = (-0.17) - (-0.20) (5.51) 

one can make a few comments: 

3 in the Standard Model 
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• the mixing with the gluon magnetic momentum operator is not much relevant. 

• The contribution from the 4 fermion operator is of the same order of magnitude, and 
with the same sign, as the one from photon magnetic momentum operator. 

• As the Ci coefficient increases with mt, the relative importance of the QCD enhancement 
decreases. 


5.4.1 Error analysis 

Let us go back to the formula in Eq. 2.15, which expresses the ratio of the b —» sy to the 
b —> ceP e decay rate. 


T (6 —> s'}) 

T (b —» ce & e ) 


6a QE D lry 

l*W’ '(!-/(Sf)) 


(mb) I 2 , 


(5.52) 


We want to briefly explain how they have been derived the theoretical error bars shown 
in Figg. 2.2, 2.5. This topic has been addressed at length by the authors of [81], and we shall 
be schematic. 

In the above mentioned papers 6 sources of uncertainty are discussed: the resulting errors 
can be added in quadrature, as there is no reason to assume they are correlated. 


(i) The use of a spectator model is now know to correspond to the leading order in the 1 /m;, 

expansion, and the existing literature [49, 65, 66] suggests that at most a 10% error can 
be ascribed to this origin. 

(ii) The cancellation of the rapidly varying factor m?, using the normalization to the semilep- 

tonic decay rate, introduces a dependence on the ratio z = m c /mb, which is known [43] 
to be z = 0.316 ± 0.013, and originates a 6% error. 

(iii) The combination of CKM parameters occurring in Eq. 2.15, assuming unitarity and 
taking into account limits from CP violation in K physics, is 


\VM\ 

IK 6 | 2 


0.95 ± 0.04 . 


(5.53) 


(iv) It is common to use a s coming from the Z physics: a s = 0.12±0.01. This determination 
however results by using NLL expressions for the observables at the Z peak, therefore 
it is in general unconsistent to use this determination at the present level of theoretical 
accuracy of the b —> sj process. The authors of [81] have checked that using two or 
one loop expressions for a s in the coefficient evolution gives an error 4 times smaller 
than the one resulting by the inaccuracy in the determination of the initial value. It 
appears therefore reasonable to use consistently a leading logarithm expression for a s 
and maintain the range of variation of the initial condition, as it comes from LEP. 

(v) In Sec. 4.2.3 we have discussed how the expression for the coefficient evolution is altered 

by going from the leading to the next-to-leading expression: we have also seen how the 
dependence on the arbitrary renormalization scale is reduced. Also in the present case it 
is possible to check [81] the explicit reduction of /./ dependence, however the central value 
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cannot be predicted without the NLL terms in the anomalous dimension matrix. On the 
other hand, as the variation of the admissible perturbative values of these parameters 
can be put in correspondence with a variation of the // scale of the LL computation, the 
resulting error can be reasonably estimated by varying the scale between, say, 1/2 mb 
and 2mb- It is worth noting that even 1/2 mb — 2.5GeV is a scale sufficiently high to 
trust perturbative QCD, so the assumption that NL coefficients remain “reasonable” is 
justified. The resulting error is about 25% 

(vi) The branching ratio BR [B —> X s ^\ is further affected by the experimental error on 
BR [B -> X c ei7 e \ = (10.4 ± 0.4) %. 

The lack of the NLL analysis is the main source of error: a NLL computation could reduce 
the uncertainty under 10%, while it is currently larger than 25%. We have already seen how 
difficult is at the moment to discriminate among the different models, in view of these large 
uncertainty. We shall come back to this point after discussing the second topic of this thesis, 
the heavy top effects at LEP, because other processes sensible to the short distance structure 
of the Standard Model, like the Z —► bb decay, put complementary limits to the New Physics 
parameters and a comparative discussion is mandatory. 




Part II 

Heavy Top effects at LEP 



Chapter 6 

Overview 


It is known from a long time that an heavy top quark results in relatively large quantum 
corrections to LEP observables. This heavy top dependence allowed to derive, together with 
low energy data, bounds on the nit value. At the present time the top quark appears to 
have been discovered at CDF, and the same calculations which helped to set the bounds are 
used to “subtract” the top effect from precision electroweak experiments, in order to uncover 
possible New Physics effects. 

In this Chapter we would like to provide an introduction to these topics and to set up a 
framework for the detailed discussion of a part of the radiative corrections, the O 
effects. 

6.1 The renormalization scheme 

It is useful to recall the methods commonly used to analyze the radiative corrections to the 
SM, limiting to the aspects more relevant for the study of the heavy top corrections at LEP. 

The Standard Model is a renormalizable theory, which means that a finite number of 
experiments is sufficient to its complete determination, allowing to make predictions for all 
the other processes. 

Given the set of “bare” parameters of the SM lagrangian, say {a*}, an equal number of 
experiments {e*} has to be performed. The theoretical relations 

e 1 (a 0 ) ( 6 . 1 ) 

can then be inverted to obtain 

<(e) (6.2) 

and specify any other experiment in terms of the basic ones. 

The choice of the defining experiments is largely arbitrary, with the only requirement of 
independence, in order to allow the inversion of the relations in Eq. 6.1. In practice, the 
accuracy in the determination of the a l D parameters is maximized choosing basic observables 
which are on one hand experimentally clean and on the other hand computable (for instance 
within perturbation theory) to an high degree of precision. 
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6.1.1 Basic observables 

The specification of the SM gauge content SU( 2) x U( 1), together with the electroweak 
symmetry breaking scale, given by the vacuum expectation value (vev) for the Higgs fields, 
require the determination of three basic quantities, which are the two gauge couplings 

9o, g'o (6-3) 

and the vacuum expectation value vq. Other than these quantities, the theory is specified 
by the fermion Yukawa couplings and the quartic Higgs self coupling. 

For what concerns the three basic parameters, an optimal set of experiments is provided 
by [17] 

• the determination of the fine structure constant a , coming from the Josephson effect or 
by the (g — 2) experiment 

a~ l = 137.0359895(61) . (6.4) 

• The Fermi constant 


= 1.16639(2) x 10 -5 GeV -2 
as determined from the lifetime of the muon and by the theoretical formula 

1 + 0 (25 - - 2 ) (l + 


(6.5) 


_i = G^rnl 
T/i “ 192?r 3 


1 - 


. m: 


mi 


a m„ 
3 7T 11 777.2 


( 6 . 6 ) 


including the QED corrections to the 4—fermion interaction. 1 
The Z— mass, as determined by the LEP experiment [90] 

M z = 91.187 ±0.007 GeV 


(6.9) 


1 One easily recognizes that the log appearing in the right hand side is the first contribution of the leading 
log series contributing to the running QED constant, as given by 




(6.7) 


where the sum runs over all the fermionic degrees of freedom which are “active” at the scale, hence one 
could have also written 


~ 192tt 3 


1 - 8 - 


1 + ±Kd ( 25 ,^ 2 ) 


( 6 . 8 ) 


which is rather an affectation at this level, but will be important when considering experiments at the scale 
of the SU( 2) x 1/(1) breaking. 
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For what concerns the Yukawa couplings, they are eliminated in favor of the fermion 
masses: we just recall that while lepton masses are unambiguously defined because these 
particles are observed in isolation , the masses of the light u. d, s quarks (or better the ratios 
m u /rrid, m s /m c i ) are extracted from pion and kaon masses using chiral symmetry, while c, b 
masses come either from the study of charmonium spectra, and by the D and B meson masses. 
In dealing with LEP observables it is possible to avoid the resulting uncertitude, for instance 
in the expression for the running QED coupling the hadronic contribution can be related to 
the e + e~ —► 7 —> hadrons using a dispersive representation. 

A different role is played by the masses of the top quark and of the Higgs field (which is a 
substitute for the quartic Higgs coupling): in the analysis of experiments up to the electroweak 
breaking scale, and in the language of effective theories, heavy fields can be removed at the 
price of renormalizing the effective lagrangian [25]: but the decoupling theorem, as we shall 
see, is “evaded”, and the corrections grow with the top mass. 


6.1.2 Derived observables 

At the tree level, the W mass is related to the basic observables by the relation 

M‘( v = M| cos 2 6w (6.10) 


where 

^e w cm*e w = - / J^. ( 6 . 11 ) 

At the quantum level, if we maintain the definition of the Weinberg angle as fixing the 
7 and Z fields as mass eigenstates, the relation with the vector boson masses is altered as 
follows 

sin 2 e w -> s 2 e = 1 - ~~^2 t ( 6 - 12 ) 

P z 

the p parameter has a tree value of 1 for a minimal SM or extensions based on doublets 
of scalar fields. However, it is also altered by quantum corrections. It is customary to define 



+ 




- s 2 
3 W 


+ c w 


A p 


(6.13) 


where the parameter s\ v = 1 — M^/M§ is just a substitute for the W mass, while Sg 
describes the coupling constants in the neutral current sector and enters neutrino scattering 
and observables at the Z peak. 

From the relation of the Charged Current coupling to the effective Fermi constant G /t 




7r a 


2 s 2 e c 2 e p M l 


(6.14) 


one obtains the definition of sg in terms of the fundamental observables, and the calculable 
p parameter 

7ra 

V2G~Mfp' 


2 2 

s e c e — 


(6.15) 
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in other words the p parameter enters the overall normalization of the Neutral Current 
vertex 




(6.16) 


and can be compactly defined as the ratio of the neutral current to charged current 
interaction at low energy 


GWc(O) 

Cecity 


(6.17) 


It must be stressed that s 2 e is not a fundamental observable, it is a derived quantity and 
can be tested in the appropriate experiments. 

The fact that the p parameter is nearly equal to one, that is, the corrections are “naturally” 
small, is the reflection of a underlying symmetry of the SM (see for instance [26]). Indeed, in 
the limit in which the JJ{ 1) gauge coupling g' goes to zero and also the Yukawa couplings go 
to zero, the SM possesses a gauged SU (2 )l symmetry and also a global SU (2 )r symmetry. 
When the Higgs field acquires a vev , this symmetry is broken down to diagonal SU (2), the so 
called “custodial” symmetry. This means that the neutral currents and the charged currents, 
being part of an SU{2) multiplet must have the same coefficients, hence p = 1 at the tree 
level. 


6.2 Heavy top effects 

The heavy top effects are a manifestation of the mass-generation mechanism in the SM, that 
is, of the Yukawa coupling to the Higgs field [10, 12], It is sometimes stated that these 
effects, growing with mt, are an “evasion” of the Applequist-Carazzone theorem [4], This 
can be misleading: simply the theorem does not apply. Recall that the theorem asserts that 
sending a dimensionful parameter, say the mass of one of the fields to infinity, the low energy 
theory results by considering the lagrangian with the heavy field removed, including finite, 
unobservable redefinition of the couplings, plus corrections which are depressed by powers of 
the large mass. But the rat —> oo limit results from a dimensionless parameter, the Yukawa 
coupling g t , going to oo. Diagrammatically, we shall see that powers of the Yukawa coupling 
from the vertices can over-compensate the m* terms from the denominators. 

6.2.1 Heavy top effects in the p parameter 

We have seen that the p parameter is protected by the “custodial” symmetry, which is broken 
by the Yukawa couplings. What really counts is not the mass of the individual quarks, but 
the mass splitting within a multiplet: for instance, if the b and t quarks were both heavy and 
degenerate, the SU (2) custodial would remain unbroken. 

The language of effective field theories is well suited to analyze the effect of the mass 
splitting [25]. Suppose we are given the t, b doublet, with mb <C M\y , mt 3> A%; we can 
divide the range of scales in three regions: below Mw the theory can be described with an 
effective Fermi theory, that is, with W, Z , t. fields integrated out. At a scale around a 
theory with the t held integrated out is used, while finally at energies of the order of mt the 
full structure of the theory is recovered. The matching between these descriptions, that is, the 
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Figure 6.1: From the effective theory to the fundamental theory 


procedure needed to derive the coefficients of the effective lagrangian, at lower energy scales, 
from the fundamental description at higher energy scales, originates the residual dependence 
of the low energy theory from the parameters of the “fundamental” theory. 

The coefficients of the four fermion operators, considered at a matching scale of about 
fj, = Mw, are connected to the gauge boson mass terms of the intermediate scale theory, while 
on his turn the gauge boson mass terms are affected by the heavy t quark present in the higher 
scale theory (Fig. 6.1). In practice then the p parameter results from the different effect of 
the heavy quark on the and the W 3 mass terms. In this language it is particularly clear 
why the effect is so relevant, and depends on the square of the top mass: the mass terms in 
the effective lagrangian of the intermediate scale theory are dimension two operators, hence 
nothing prevents them to depend on m|. 

We anticipate the result at one-loop order, 

' =1+n '£$ (6A8) 

where N c is the number of colors, which shows that to this order the effect comes entirely 
from fermion loops. 

6.2.2 Heavy top effects in vertices 

The effect on the p parameter can be regarded as a self-energy effect; in the scheme which 
sets among the fundamental parameters the Z mass and the Fermi coupling G it results 
from the W ^ self-energy. 

Other heavy top effects stem from the vertex corrections, and the most relevant, from a 
numerical point of view, are the corrections to the Z —> // vertex [33, 39, 40, 41]. We shall 
see how in the limit of a very heavy top only the exchange of the would-be-Goldstone fields 
is relevant to the computation of the corrections, as in Fig. 6.2. 

The vertex is proportional to |/1 2 , hence due to the smallness of Vtd, Vts coefficients 2 
the heavy top effects are relevant numerically only for the Z bb vertex. In the following we 
shall always set Vp, = 1. 

Due to the SU (2 )l symmetry, a purely left handed correction results to the lowest order 
Z —> bb vertex: we again anticipate the result for the amplitude, at leading order 


2 \V t d\ = 0.004 A 0.015, \V ts \ = 0.030 4-0.048(70] 
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Figure 6.2: Leading corrections to Z —► bb in the mt —» oo limit. 
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(6.19) 


Apart the overall factor, the corrections can be included in a redefinition of the vector 
and axial electroweak couplings, 

4 9 , 

g w = 1 - -sg + r g bA = 1 + r ( 6 . 20 ) 

thus separating the vertex corrections to the amplitude, included in r, from the corrections 
induced by the very connection to the fundamental observables, included in s$. 


6.2.3 Bounds on top mass 

We have to remind, before discussing the bounds, that a comparison with experiment is 
possible only introducing the main sources of radiative corrections to electroweak observables, 
coming from the “pure” QED effects 3 and by the QCD corrections. 

Their inclusion is mandatory and straightforward: in this thesis it is impossible to account 
for the whole set of radiative corrections, but it goes without saying that in graphs and tables 
all the relevant effects other than the heavy top ones are included whenever needed. 


Low energy data 

The most precise low energy informations come from the determination of NC/CC neutrino- 
nucleon cross section ratios [ 68 ], by the CDHS[35, 45] and CHARM[34] experiments: the 
accuracy is about 1 %. At tree level, this ratios are given by [23] 


(M w y l-2s 2 w + f (1 + r) si 
2(l-s 2 w f 

_ fM w \ 4 1 “ 2s w + T ( X + y) s w 
VMzJ 2(1-S 2 W ) 2 


( 6 . 21 ) 


3 As shown by Okun and collaborators [91, 108], for some time QED effects have been sufficient to account 
for all the radiative corrections. 
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where r = accl G CC — 0.4. Leaving aside the QED, box and vertex corrections, the large 
mt effects can be incorporated in this formulas by replacing s'f v —> s'f v + A pcfy. But in the 
expression for R u the dependence on s/ v is weak: 

Rv = 9 ((i + °-56svy + O • ( 6 . 22 ) 

Hence the ratio Mw/Mz, when extracted from R u , depends very weakly on rrq! this 
stability is the result of an accidental cancellation, and does not appear neither in the Rp 
ratio nor in the v e cross sections. Combining this Mw/Mz measurement and the other 
low energy u e, u q, eq data one can obtain an upper limit nit < 168GeV at 68 %C.L., for 
Mu ~ Mz[ 51], but including also the direct measurement at CDF of the ratio of vector-boson 
masses, together with the theoretical prediction from the input parameters a, GMz , one 
can obtain a global fit (see for instance [68]) 

m t = 122^3® GeV (68%C.L.) (6.23) 

for 40GeV < M„ < ITeV. 

LEP data 

We refer to a number of good reviews [44, 53, 60, 67, 92, 95] for a complete discussion of 
the observables studied at LEP around the Z resonance. We shall here be brief because we 
will come back to this subject in the Chapter 8, when also the two-loop corrections will be 
considered. 

It suffices here for the purpose of illustration, to focus on one of the most precisely mea¬ 
sured observables, the total width T z- In Fig. 6.3 it is shown the comparison between theory 
and experiment. The various curves, plotted in dependence of mt, for different values of the 
Higgs mass, are superimposed on the band of experimental values. 

One can note a few features common also to similar plots for the various partial widths, 
like T (Z —► p + p~) or the forward-backward asymmetries: the Mu dependence is rather slow, 
as a consequence of the Veltman “screening” theorem [9]: leading quantum corrections depend 
logarithmically on Mu- 

Combining limits extracted from plots like the one in Fig. 6.3 with low energy informations, 
and the determination of M\y from CDF, different [68] analyses agree on the following limits 
on mt, at la level 

m t = 144l|l^GeV (6.24) 

where the second error is desumed by varying Mu around a central value of 300GeV 4 . 

The very significance of these bounds relies on the stability of the perturbative expansion: 
predicting an heavy top mass starting from an expansion in powers of mt calls for a check of the 
next perturbative order. In the following we shall describe the calculation of gf effects to LEP 
parameters, which has been useful, at the time it was performed, to improve the confidence 
in the existing results. Presently the gf effects are included in the analyses of experimental 

4 A more recent analysis [111] finds a result in remarkable agreement with the recent CDF results, namely 

mt = 175 ± llt^GeV for M H = 300t^GeV . 
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Figure 6.3: The total Z width from LEP 
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data, since some observables are being measured at a comparable level of accuracy: however 
we shall see that, in view of the present probable value of the top mass, the effect cannot 
be considered dominant, and eventually a full two loop calculation shall be mandatory. 








Chapter 7 

Strategy to compute heavy top 
corrections 


Our aim is to evaluate all the corrections growing like G lt rrif , G'j t rnf to the Standard Model 
parameters 1 ; we have seen that it is known how these effects stem from the Yukawa sector of 
the theory [52, 50]. The idea is to define an appropriate limit such that for g t 2> g, g' the theory 
reduces to a pure Yukawa model, that is, to a scalar theory, simplifying the computation. 

The easiest way to define this limit is to work in a renormalizable gauge: in this way only 
charged and pseudoscalar Higgs fields, not “eaten up” by the gauge fields, have couplings with 
the (i, b ) doublet proportional to gt- As a consequence only Feynman diagrams resulting by 
the exchange of these scalar particles are to be considered, while the quantum gauge fields 
are irrelevant, and the gauge fields appear only on the external lines and can be considered 
of classical nature. 

These considerations are made rigorous by using the background field formulation of the 
Standard Model, which allows to relate different Green functions by classical Ward Identities; 
in particular it results that quantum corrections to the p parameter and to the Z —► bb vertex 
are connected to correlations of scalar fields. The advantage over a standard formulation, 
which would lead to Slavnov-Taylor Identities, is that it is made manifest how these large gt 
effects are completely decoupled from the gauge structure of the theory. The gauge structure 
almost completely disappears, in the sense that gauge fields do not propagate, and they couple 
as external sources to the classical SU( 2) <g) U( 1) currents of the theory. 

It can be useful to give a pictorial resume of the strategy 

Full Model 

Classical SU (2) ® U (1) gauge symmetry —> BRS symmetry 


Fields 


(A, Z, W) U (H, X , U (t, b) 


Inputs 


&e.m. G> M z m (m H ) 


Mimes functions of the ratio mt/rriH, and disregarding QCD effects 
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4 

Background Field Gauge 

• Fields doubling <h —► + <h el 

• Ward Identities on external lines 

• Slavnov-Taylor Identities, ghosts, ... in internal lines 


4 

9<-9t —» M W) z < rri t 

ghosts and gauge fixing related vertices irrelevant; gauge bosons coupled as external currents 


t 


Gaugeless limit of the Standard Model 


7.1 The Ward Identities 

In App. C and D the formulation of the quantum Standard Model in the renormalizable 
non-linear gauge and then in the background field gauge is discussed. 

The formulas given there are rather involved, in the attempt to write down results as 
general as possible, but the philosophy is quite elementary, and the results are self explanatory. 

Having chosen to preserve the classical symmetries of the theory, by using a background 
field formulation, we are allowed to impose on the renormalized effective action T the same 
Ward identities 2 we would obtain in the “tree-level” theory (see Appendix D). 

As discussed in the previous chapter, for the purpose of evaluating the large rrq effects we 
can reduce the matter content of the Standard Model only to the third quark doublet ( t , b). 

By exploiting the invariance under the so called type II transformations, which are pure 
rotations on the quantum fields, we can obtain 4 generating equations for the classical Ward 
Identitities relating 1 PI Green functions, corresponding to the 4 generators of SU( 2) <8> U[l) 

E 5 Vf Cl (*) JJ, = 0 * = 1 • • • 4 • ( 7 - 1 ) 

Taking into account the formulas C.39, 0.40, C.41 we can easily write down their detailed 
form: in the following the =F, ± symbols whould be understood as summations over the 
contributions from oppositely charged fields. 


A-current Ward Identity 

ST 


d x - 

SAfj (x) 


, ST , <jr 

= ±ieW± _ l ± ie^ 


M 6W± 


5(j) d 


^ 1 idQq 

q=t,b 


ST 

5q(x 


■q(x) + q{x) 


ST - 
Sq(x)_ 


2 in absence of external quantum legs or ghosts 


(7.2) 
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Z-current Ward Identity 


d: 


<5r 


^SZ^x) 

Pig cosri/ W, 


ST 


le 




ST 


+ 


M 8W± tan ‘20ii/ S<iP 1 sin 20 w A SH 


e c>r e 

XTtj -( H + V )Z; 


5T 


-E- 


le 


smvF cosyp 


ST 

Sq(x) 


sin 2 9 W Sx 

ST 


(T 3 P L - sin^/ Q) q(x) + q{x) (t 3 Pr - sinjv Q) 


5q{x)\ 

(7.3) 


here the projector Pr on left states has been explicitly introduced. 


Charged current Ward Identities 

Here no summation on signs ±, =F is understood 


(d,, ± igW ; 


8T 

~5wZ 


^(e^+gf 


SA„ 


SZ„ 


g,= ( . ST 5T 

+ ^ f * 1 sh + K 


+ 2 (±HH + v) -x) 


5T 

JZ± 


w \' ( ± dr \ 

— - : y I - rr ctT 4- nr rr - I 


a/2 q \Sq L 


<y 9 l + QL& 


^qlJ 


(7.4) 


Having fully described the symmetry content of the theory, we can specialize to the pro¬ 
cesses under investigation. 


7.2 Identities relevant for Z —> bb and the p parameter 

By differentiating Eq. 7.3 with respect to b fields one easily obtains 


d\ 


ST 


** 8Zp(x)b{y)b(z) 


ev 


5T 


sin 2 Ow 8x{x)5b(y)Sb(z ) 


*e 


Slllw cosvf 


8T 


5b{x)5b(z) 


~\ Pl + 3 Shl ^ ) 6 ( X ~ y ^ + 


+8{x - z) ( ~Pr + ^ sin^ 


ST 


5b(y)8b(x)\ 


(7.5) 


which allows to express the irreducible Z —> bb vertex in terms of the corresponding vertex 
for y —> bb 3 . 

3 Note a useful identity relating the W + —> tb vertex with the corresponding one for (j> + 
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The identities relevant for the p parameter can be analogously obtained by considering 
the wave function renormalization of fields Z, x and W±, 4>^ 4 and the tadpole term for H: 


d ^8Z^x)6Z u (y) 
ox W 

»6Z^x)5 x (y) 


gv ST 
2 cSx(x)5Z u (y) 

9 \ xf ST 1 

2 c[ X V 5H(x) V 5x(x)Sx(y)_ 


(7.7) 


As a renormalization condition we shall set to zero the tadpole term, thus imposing that 
the vacuum expectation value v does not renormalize: consequently the two relations in 
Eq. 7.7 combine to give 


QXQV 6T = 9 2 V 2 8T 

** lt 8Z li (x)8Z v {y) 4 C 2 8 X (x)8 X (y) ' 

We can also relate the renormalization of W ±, (j)^ fields by writing 


(7.8) 


d ; 


ST 


d\ 


M SW+{x)S(t>-{y ) 

ST 

»8W+(x)SWZ{y) 


= o 8(x-y)[-i 


= iS(x - y) 


ST 


+ 


ST 


SH(x) Sx(x) 

ST ST 

e „ . , ' + g- 


+ 


+ 


igv 


8T 


2 8<f>+(x)<f>~(y) 
igv <5T 


L 8A v {y) 8Z u (y) \ 2 6<J>+(x)8Wv(y) 


(7.9) 


and then, again setting to zero the tadpoles, 


d x d y _—_ 

M u SW+(x)SWu(y) 


gv_ 

2 


ST 


5<f>+(x)(f>-(y ) 


(7.10) 


7.2.1 Form factors 

In momentum space, the identities in Eqs. 7.8, 7.10 can be rewritten as follows: 


PnVvTi z ^{p) = n x (p) 

PuPutfOib) = (jpj n <p(p) . 

The self energy at zero momentum gives the correction to the mass, 


(7.11) 


ar _ 

SWf (x)8b(y)Si(z) 
_ig_ \ ST 
72 St(x)St(z) 


igv <5r 

2 5(j>+ (x)Sb(y)5t(z) 

P L S{x -y) + 5(x - z)P R 


8T 

8b(y)5b(x) 


(7.6) 


4 Note that the background gauge fixing sets to zero the mixing Z , \ only in the quantum quadratic part of 
the action 
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n %, = (-*) VnuSm 2 z 

p =o 

= im ( 7 . 12 ) 

where the mass mz appearing in the formula is the tree level one. 

We shall use the notation, for the self energy of the scalar fields, 

11/ = ip 2 S/ + O (1) 

and obtain the mass renormalizations 


S m 2 z 

,dli x 

m 2 z 

1 dp 2 p=0 

Smlv 

.oTF'l 

m w 

1 

to 

43 

II 

O 


The unrenormalized ratio p is then given by 


(7.13) 


_ _ rrtfy + 5 rrtfy 
^ c (m 2 z + Srri z ) 

= i +«(4" - 42) + fi2 [(4 2) - 4*) + 4 1 ’ (4" - 42)] < 7 - 14 ) 

where the orders of Ti have been put in evidence. We shall see in the following that a term 
(l + S^±) will be readsorbed by the renormalization procedure in the definition of the Fermi 
coupling G p . Consequently the renormalized p parameter will be obtained by the formula 

p = i+* (4" - 42) + fi2 [(4" - 42) 2 + (4 2) - 42)] + 0 (* 3 ) • < 7 ' 15 ) 

which deserves a few comments: first of all, it is made evident that only the difference of the 
self energies for charged and pseudoscalar unphysical Higgs fields has to be computed. Next 
we note that at h 2 order there are two contributions, one resulting from reducible diagrams 
and proportional to the square of the one-loop coontributions, the other coming from two-loop 
diagrams. It is therefore common to present the results of a perturbative evaluation of the p 
parameter in a different form, 

i - 1 = -ft (EW - E«) - ft 2 (E« - EW) + O (ft 3 ) . (7.16) 

A similar analysis for the Ward Identity in Eq. 7.5 (relevant for the Z —> bb vertex) is 
slightly simplified if we choose as renormalization condition for the b propagator a unit residue 
at the mass pole (at mb = 0 in the present approximation): the tree level terms cancel in 
both sides of the equation, and it is immediate to show that the quantum corrections of the 
vertices at zero exchanged momentum are related by 

( Z ^) quantum = "^7^ ||fc=0 (7-17) 

where k is the momentum entering the x> Z line. 
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7.3 The Yukawa model 

The computation of the p parameter and the Z —» bb vertex have been reduced, through 
Eq. 7.15 and 7.17, to the evaluation of correlations of scalar fields. Now it comes the final step, 
in the heavy top limit these fields are described by the following Yukawa lagrangian, resulting 
by dropping the gauge couplings and fields from the SM lagrangian, in the renormalizable, 
background field gauge fixed form: 


-C = —Ql<PQl ~ tnfttn - bfifthR - ( d(d M $) 

- V (V<h) - g t {Q L *t R + i&Q L ) 

Ql = (l) * (718) 

with the usual form for the potential V 

V (. x) = p 2 x + Xx 2 . (7.19) 

The Ward Identities we have discussed before are now derived by supplementing this 
lagrangian with the coupling of the external W, Z fields to the conserved classical currents, 

A£ = ~T G 7 -" 2 7 m ) - d (Y K + w) • (r.20) 

resulting from the different symmetry generators 

• baryon 17(1) 

Jfl = V + bub (7.21) 

• weak hypercharge 17(1) 

J} = b’L + T^rr/iJR + ^>LhJ>L - ‘^bra/ibn - $*8^® + (7.22) 

• T 3 generator of 517(2) 

Jl = ~ hlbL + b + d^~ - ft-d^ + ) + (<^$>0 - <M^S) (7.23) 

• Combining in the usual way the currents J 3 and J 5 one gets the e.m. current J em 

J e b = ~ \b/* b + b+dub - ft~d fl ft + ) (7.24) 

• T± generators of SU (2) 


= *l7/A + -j= -ft dpH) + -^= (xd^ft -ft d^x) 

7+ = -bLX^L + -^= (Hd /1 ft + - ft + 8pH) - -b (xdnft + ~ ft + 8 fi x) • (7.25) 
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In other words, the classical Ward Identities reduce to the conservation of these Noether 
currents; in particular, as the electromagnetic current is exactly conserved, we have that only 
the renormalization of current J 3 is important. 

Note that Eq. 7.18 is written without specifying the phase of the theory (this is the 
meaning of the H' notation): when considering the symmetry breaking a subtlety arises and 
the conclusions of the previous section, that is, the Ward Identities relating the effective Z 
and x vertices, are to be translated in this language by stating that the current coupled to 
the Z field is still 

J* = tLhxtL - bLhxbL + (<f> + d,j,<l>~ - 0 - <9 m 0 + ) + i (Hdf.x ~ xd» H ) (7-26) 

but the conserved current is 

4 = jfi+ivd^x, (7.27) 

meaning that, as operators 

dfj Jn = — iv^x + motion equations (7.28) 

and the insertion of the current J^ is put in relation with the insertion of operator d^x- 
Analogously one has that 

4 = 4 + (7-29) 

are the charged conserved currents. 

We stress that this is an information coming from the full model, in the non-linear back¬ 
ground field formulation, which allowed to derive the Ward Identities. 

In summary, by considering the SM quantized in the background field gauge, we have 
shown that the same Ward Identities, in the complete theory and in the so-called gaugeless 
theory, relate the quantities most sensitive to mt to correlations of scalar fields, which can be 
studied in a pure Yukawa model. 




Chapter 8 

Heavy Top effects in the Gaugeless 
Limit 


We are now in position to describe in detail the computations needed to obtain the O [G^mf j 
corrections. We start by defining the renormalized model. 

8.1 Renormalized model 

We work in dimensional regularization using the NDR scheme for the 75 matrix. We have 
seen in considering the computation of QCD corrections to the b —> sy process that up to two 
loops it is possible to avoid any resulting regularization ambiguity. 

The following renormalization constants are introduced, to cancel poles and perform ad¬ 
ditional subtractions in order to impose suitable conditions, without breaking the classical 
symmetries 


Ql 

- z' l /2 Q l 

t R - 

- 4 /2 *« 

at - 

> Zg9t 

$ - 

► z\ ,2 <$> 


h 2 - (m 2 - 5 ^) V 
A - A Z X Z^ 2 
v —> z]J 2 (v — 6v) 


( 8 . 1 ) 


where the redefinitions of v, Z\, /i 2 contain also Z^ for reasons of opportunity. Indeed, one 

has 


( \ 7 l /2 ( ^ \_ 7 l/ 2 . 

I H+ix _) _ j ^<t> l H +*X _|_ v-Sv J — 9 . 

\ \/2 \[2 ) \ \[2 \/2 / 

The renormalized lagrangian for the Yukawa model is therefore 


( 8 . 2 ) 


C = -Z L Q L pQ L -Z R t R pt R -b R pb R -Z 4 > {d l &) li {d l &) 

- g t Z g {Z L Z R Z (j> ) ll 2 {Q L ^t R + i^Q L } 

- (V - 5/u 2 ) $ - AZ a (V<h) 2 


(8.3) 
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where for each Z constant we mean Z = 1 + hSZ. No wave function renormalization is 
needed for the which completely decouples in this reduced model. 

Apart the subtraction of poles, the counterterms arising from these redefinitions are set 
according to the following renormalization conditions: 


Tadpole cancellation 


Higgs mass at its classical value run = 2Xv 2 


T + iv 3 ^ r + 2A— -XZ x =0. 


I p 2 =—m 2 H 


+ i dfj 2 — 3v 2 X 5Z\ + 6A v5v = 0 


Higgs Propagator residue set to 1. In order to preserve the explicit invariance of the 
action only one of the residues of the Higgs fields can be fixed; choosing the physical 
one, the condition is 


d ^ 

o 9 


— idZrh = 0 . 


OP I p 2 =—m 2 H 

Masses of Xi 4^ ■ Hi this restricted model we forget masses proportional to the gauge 
couplings, when compared to the nit, run, so we set to zero the masses of the unphysical 
degrees of freedom. 

+ i 5[i 2 — v 2 X5Z\ + 2Xv5v = 0 ; (8-7) 

p 2 =o L J 

it is crucial that the two self energies are equal at p 2 = 0. As one may expect this 
condition coalesces with the one on the tadpole and no new constraint arises. 

Bottom propagator . In this model it is easy to see that no mass is generated for the b 
quark, so the only condition we impose is on the residue 

n 6 = 

Xj; + 6Z L = 0 ( 8 . 8 ) 

Top Mass . The perturbative expansion generates a non zero axial term for the top self 
energy, which can be set to zero at the physical pole: given the expression 


n* = 75 + im t T, s 


one can impose 


+ - ( 5Z l - SZ R ) = 0 


(8.10) 
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We fix the mass at the tree level value with the condition 

\i> 'Ey + — (5Zi + 5Zr) + imt Eg -> = 0 ) ( 8 - 11 ) 

l L 2 j L m J J j 6 =imt 

where Suit is only a shortcut for 

1 s 

5mt = 6Zg + — (5Zl + 5Zr + SZ^) -; ( 8 - 12 ) 

it follows that 

( 1 \ s 

6Z g + —5Z,pj H-= 0 . (8.13) 

L z J V 

Decay constant G^ . The value of the fj, decay constant cannot be of course obtained by 
this model: we may use an indirect argument, noting that contributions to the process 
/i —> e u e Ug arise from corrections on the line of the (fZ exchanged: these corrections, 
together with the mass and wave-functions renormalizations, give rise to a non zero 
residue which can be incorporated in the definition of the physical coupling. 

In short, if 

n0± (p ^ — n 0 ± (0) = ip E 0± (8-14) 

is the (f>^ self energy after mass renormalization, then the propagator has the residue 




again terms proportional to the gauge coupling g have been omitted. 
The resulting correction on Gg is given by 


(8.15) 


G l loop 




Now the tree level relation between Gg and the other couplings is 


(8.16) 


G± = J_ 

y/2 2v 2 


(8.17) 


so we can proceed by fixing also the residue of the <$£ propagator to 1 and setting to zero the 
correction to v, thus preserving the relation 8.17. 

It is to be stressed that the only physical inputs are nit , mu- The other conditions 
like the ones on the propagator residues are just conventional, to simplify the expression of 
amplitudes. 
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Figure 8.1: Legenda for fields in Feynman Diagrams 


H 



I 

I 
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Figure 8.2: Higgs “tadpole” diagrams 


8.1.1 The renormalization constants 

Here we list the one-loop computations needed to set up the renormalization framework up 
to two-loops. 

The abbreviation 

A =-7£ + log47r (8.18) 

is used, where e = (4 — d) /2, and the fields are addressed according to the legenda in 
Fig. 8.1. 

Tadpole (Fig. 8.2) 


T 

v 


3 m 


H 


(47r) 


A + 1 — log 


m 


H 




— 2m 


2 9t 


(4tt 


~2 N c 


A + 1 — log 


m 7 




where N c is the number of colors. 


(8.19) 



Figure 8.3: Higgs or x self energy 
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8.2 One loop order 
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H,X 



Figure 8.6: Top self energy 


t self energy (Fig. 8.6) 


n t = S vi> + 75 + irnt^s , 


'V 


Sa 


p 2 =—m 2 


p 2 =—m 2 


p 2 =—m 2 


St 
(47r) 

St 

(4tt) 

St 

(4tt) 


3 11 1 3 mf 

I A + T + 4 l0Er -4 l0S 7? 


1 . 1 1. m? 

-4 A -2 + 4 l0g 7? 


1 1, to: 

-«+ o 1Q g 


m 


The renormalization constants are then 


6Z l = 
5Z r = 

6mt 

rn t 


St 


(47r) 

St 
(47r) 

(4t r) 


3 1 . 1 n TO: 

-A + - log — 

4 2 9 b fi‘ 


H 


7 


9' 

— j - A + i°g —f- 

4 ju 2 


3 3 

-8 + 4 l0gr 


(8.23) 


(8.24) 


The 4i; constant is set equal to zero: it follows that the counterterm Srrit/mt defines also 
the finite vertex counterterms, in accordance with Eq. 8.3 and Eq. 8.12. 


8.2 One loop order 

The expression for the (j)±, x self energies in Eq. 8.21 together with Eq. 7.15 allows to obtain 
the p parameter at one loop order 


p=l + 5p\ = l + 


N C St 


2 (47tF 

where N r = 3 is the number of colors. The tree level relations 


(8.25) 


St = 


\[2 mt G 


M _ 


1 


y/2 2 V 2 


(8.26) 
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I 



Figure 8.7: One loop xbb vertex 
can be used to write the result in terms of the physical inputs, 

'= i+ "'£$- (8 - 27) 

The computation of the single Feynman diagram contributing to the xbb vertex, in Fig. 8.7 
allows to obtain 

= (8.28) 

and using Eq. 7.17 one obtains the quantum contribution to the amplitude 

quantum = ^ 

Before expressing also this result in terms of physical inputs, let us consider the two-loop 
contributions. 


8.3 Two-loop order: Z —> bb 


The Feynman diagrams contributing to the (xbb) are listed in Fig. 8.8. The results are 
normalized to the one loop calculation, and an overall factor 


90 = ~ 


H L 9t 
v (47r) 2 


(8.30) 


is understood, where q is the momentum entering the x vertex. 

Only two mass scales are present, m* and m#, and apart an overall scale factor and the 
renormalization logarithms they appear only in the combination 


mi 


r = 


TYl 


H 


(8.31) 


The “technology” for the computation of this diagrams is well known, and we have ob¬ 
tained analytical results for arbitrary values of the ratio r; for compactness, the results graph 
by graph will be reported in the asymptotic small r limit, and when appropriate in the large 
r limit, but the full analytical results will be also discussed in the following. 
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Two loop contributions to the vertex 


Non divergent diagrams 

We recall that the results are normalized to go: 


Diagrams 1 — 4, 7 


9i = 


92 = 


9t 


(47r) 
9t 
(47r) 


1 3 , 

-1— log r 

16 8 5 


r N r . 
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94 

97 


9t 

'27 

7T 2 

11 

( 47 r ) 2 

32 

12 

16 

9? 

1 

7r 2 

7 

( 47 r ) 2 

16 

+ 12 

8 

9t 

"l 

7r 2 

Ilof 

( 47 r ) 2 

[ 4 + 

24 + 


log r -|-log 2 r + O (r log r) 

16 


(8.32) 


Divergent diagrams 

It is convenient to impose some of the renormalization conditions diagram by diagram: in 
particular, the mass renormalization of the (fr^ fields entering the diagram 5. We list also the 
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counterterm for the t self energy entering graph 6 . 


Diagram 5 The result after mass subtraction only is 


^mass subtracted 

9 5 


Y - ^ - p^logr + O (rlogr) + 


+ - 


(4 ttY (47r) 




(8.33) 


where as usual A = ^ — 7 e + log 47 r. The wave function renormalization 5Z^ coming from 
the condition on the Higgs propagator residue does not cancel the second line in Eq. 8.33 
completely: it leaves a finite result which will cancel against the redefinition of G^. 

Diagram 6 Here we show separately the graph and the counterterm induced by mass and 
wave function renormalization. 


^ [4 4 +^ + £++ 2iog ^ • 

The counterterm results from the insertion of the following self energy correction 


(8.34) 


\ 3 1 . 1 . to? 

Ct = * -4-2 A+ 2 log 7? 


+ i> R -- — A + log ■ 


"3 3 

+ im t - - - log r ; 

when inserted in the one loop graph the counterterm results 


1 1 Q Q 

counterterm = Jh _ _ 3 + * log r - 2 iog 

(47r) [ e 8 4 M 


The total result is therefore 
gt \ 9 7 r 2 1 

96= (iW[8 + 24-2 l0gr _ ' 


(8.35) 


(8.36) 


(8.37) 


Other finite counterterms 

The counterterms for the t and self energies have been sufficient to obtain a finite result: 
it only remains to include finite counterterms resulting from the Yukawa couplings. Having 
fixed Sv to zero, as a consequence of the definition of G the counterterms for the vertices 
(xtt) , (0 + 6i), • • • are connected to the mass counterterm of the t: each of these vertices is 
multiplied by the factor 


c - 

1 + 5Z g + — (5Zl + 5Zr + 6Zj) = 1 + ~y~ > 


(8.38) 
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hence one has a contribution, proportional to the one loop result, of the form 
,5m t 


9c = 3- 


m t 


9t 

(4vr) 


9 9 

-g + 4 lo g(0 


(8.39) 


8.3.1 Cumulative result 

It is convenient to define the dimensionless expansion parameter 

x = G n™t 
8tt 2 \/2 

and write the (x bb) in the form 

{X»> = *r 

V 

where 


(8.40) 


(8.41) 


r = — 2x (l + . 

In the small r, (m# 3> rrit) limit we have 

_C 2 ) 311 + 247T 2 + 282 log r + 90 log 2 r 

1 = 144 


(8.42) 


(8.43) 


while in the large r limit 

r (2) = . (8.44) 

To understand the meaning of the r parameter, note that the Z —► bb amplitude can be 
written as follows 


41 (Z„ - bb) 


9v 

9a 


4 c 9a r y^'y§) 

4 9 

+T 

1 + r . 


(8.45) 


We recall that s = sin 6\y is not in our renormalization scheme a physical input and it is 
eliminated through the relation 1 in Eq. 6.11 


% = \ ( n / 2 G„MIp) ■ 

V should not be confused with the definition of 


sin^ 


M 2 Z 


(8.47) 


(8.46) 
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The decay width is therefore given by 

r (Z - bb) = ^0Mlp ( \g v \ 2 + \g a \ 2 ) , (8.48) 

where contributions proportional to mb/Mw have been dropped. They will be reinserted 
when dealing with numerical results and physical observables. Here it was interesting to 
show the dependence on the p parameter, coming from the renormalization conditions. This 
dependence is conventionally called oblique , as it is common to all the physical observable as 
a consequence of the theoretical framework, in opposition to the non oblique corrections, like 
r itself. 

We require now the full two-loop heavy t expression for the p parameter. 


8.4 Two loop order: p 


Let us recall that defining the renormalized model in Sec. 7.2.1 we have fixed G^ by adsorbing 
in the physical value the self energy contribution from c f>±: we have used the relation 


^yphys _ ^fbare 


b+ 

Explicitly, in the formula in Eq. 7.14 


^ = (1+E<2)(e?> 


y (4) 

- 


+ 


( 4 1 * 




+ 


( 4 2) 


- E 


(8.49) 


(8.50) 


the prefactor in the first term was adsorbed in the definition of Gp. and bp results calcu¬ 
lable in terms of finite measurable quantities 2 . 

The two loop graphs needed to compute the difference of self energies of the unphysical 
Higgs fields at two loop level are listed in Fig. 8.9, and the individual contributions are 
addressed as follows 


bp# = E x - E^ (8.51) 

where # numbers the graph, starting from 2 because the number 1 was used for the one- 
loop result. Again the detailed results are given in the small r limit, but the full analytical 
formulas have been used in the numerical evaluation. 


Diagram 2 There is no contribution to (j) se lf energy, while the x self energy gives 


bp2 



\ + ~F + \ ( logr + log2 r ) 


(8.52) 


2 We include in the definition of Sy the subtraction of counterterms, both infinite and finite parts, as in the 
case of the Z —► bb process 
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Figure 8.9: Two loop graphs for the computation of p parameter 

Diagram 3 In the difference of the self energies the square pole cancels 
Spf Te = 
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(47r 
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The t self energy counterterm gives the contribution 
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(An 


2 N c 


7 3/1 , ml 

8 + 2(Te-' ,E - loe l^ 


and the sum is finite as expected 


(8.53) 


(8.54) 


5/03 = 


9t 


,(4tt) 


N r , 


39 7T 2 1 1,2 

-1— log r -log r 

16 4 8 s 8 s 


(8.55) 


Diagram 4 

dpi = 


Jt \ [3 7T 2 

,(JT j c (2 + T + logr y 


2 


(8.56) 
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Diagram 5 


$P5 




(8.57) 


Finite counterterms As in the case of the Z bb process, we have to include finite vertex 
counterterms, in analogy with Eq. 8.39 


$P6 = 


25m/ 


m t 


Spi 


= 2 


9t 


+ 7 lo g r ) 6 Pi 

2 


(47t) 2 V 8 4 


9t 


3 3 

-g+ i 1 °gr 


V( 47r f j 

The two loop contribution is therefore 


(8.58) 


4 21 - = E S »r, 


n =2 

g 2 V at 4 ^ + 471-2 + 54 log r + 6 log 2 r 
2 


\( 47r 

We already know that 

SW - Eg = Sp i = —, 

x ° ' 2 (4vr) 2 

and combining the results we obtain 


16 


(8.59) 


(8.60) 


5p = hSpi + h 2 

N c 9t 


5p\ + 5p,, 


n=2 


2 (4 


:7T 


IT 


9t_ 

8 (47r) 


^49 + AN C + 47t 2 + 54 log r + 6 log 2 r^ 


= Af r .x 


1 + x 


49 + 41V C + 47r 2 + 54 log r + 6 log 2 r 


(8.61) 


This result confirms a preceding computation [36], performed in the mu = 0 limit. 
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8.5 Summary of the results 

We have already discussed in the previous chapter that it is common to present the p param¬ 
eter as 

i-l = -JVc*(lW 2) ) < 8 - 62 > 

recalling the definition of t 

t = — 2x (l + , (8.63) 

the asymptotic expansions of the two loop contributions p ^, r® are 

(i) for r = ( mt/mii ) 2 1 

(2) 49 2 27 , 3,2 

p ( ’ = — + TT + — log r + - log^ r + 

+- ^2 — 12-7T 2 + 12 log r — 27 log 2 r'j + 

2 

(l613 - 240 t r 2 - 1500 log r — 720 log 2 r'j + O (r 3 ) 

(8.64) 

t ( ' 2 ^ = 311 + 247T 2 + 282 log r + 90 log 2 r+ 

—4r (40 + 67 t 2 + 15 log r + 18 log 2 rj + 

(24209 - 6000tt 2 - 45420 log r - 18000 log 2 r'j + O (r 3 ) 

where some sub-asymptotic contributions are included in order to improve the convergence 
to the exact formula. 

(ii) for r 3> 1 

p« = 19-2^-11+ 0(1211) 

v r V r ) 

(8.65) 

r» = 27-^ + 4^/logp, 

3 yjr \ r ) 

The r —> 00 limit of eq. (8.66) confirms the result given in ref. [36]. 

These asymptotic expressions are not sufficient for an accurate evaluation of the second- 
order coefficients p^ and r*- 2 ^ for mt — mu- 

In Tab. 8.5 the numerical values of p^ and r ^ as functions of rrin/rrif are therefore 
given. 

The exact numerical results are compared with the asymptotic expansions in Fig. 8.10: 
one can see that already for mu > 2 mt the asymptotic expansion is fairly accurate. 
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m H /m t 

pGi 

tV) 

0.1 

-1.82 

4.67 

0.2 

-2.70 

3.90 

0.3 

-3.46 

3.30 

0.4 

-4.13 

2.83 

0.5 

-4.72 

2.46 

0.6 

-5.25 

2.16 

0.7 

-5.74 

1.92 

0.8 

-6.18 

1.73 

0.9 

-6.58 

1.59 

1.0 

-6.95 

1.47 


m H /m t 

pGi 

tV) 

1.1 

-7.30 

1.38 

1.2 

-7.62 

1.32 

1.3 

-7.91 

1.27 

1.4 

-8.19 

1.24 

1.5 

-8.44 

1.23 

1.6 

-8.68 

1.23 

1.7 

-8.90 

1.24 

1.8 

-9.11 

1.26 

1.9 

-9.30 

1.29 

2.0 

-9.48 

1.33 


m H /m t 

pO) 

tV) 

2.1 

-9.66 

1.37 

2.2 

-9.81 

1.42 

2.3 

-9.96 

1.47 

2.4 

-10.1 

1.53 

2.5 

-10.2 

1.60 

2.6 

-10.4 

1.66 

2.7 

-10.5 

1.73 

2.8 

-10.6 

1.80 

2.9 

-10.7 

1.87 

3.0 

-10.8 

1.95 


Table 8.1: 


Values of pG) and rG) as functions of rn ^/rn t 


These results have been confirmed by other authors [96, 106]; they have also been able 
to fully exploit algebraic identities of the Li 2 functions in order to obtain compact analytical 
results 3 . 

All in all, the corrections to the one-loop computation are small, a fact that justifies a 
posteriori the perturbative expansion, even for rather large values of the gt coupling. Indeed, 
for mt — 180GeV one has gt — 1.0, and the typical expansion parameter 


9t 

at = 4 tt = 


GiJn't 


7r 


V2 


( 8 . 66 ) 


has the value of at — 0.08. 

The good convergence of the perturbative expansion is shown in Fig. (8.11), where the 
first and second order expressions for 4 — 1 and r are compared, for different values of mu 
in the mt = 150-250 GeV interval: the most affected quantity is the p parameter, and even 
for mt = 250 GeV one has less than a 10% correction of the one-loop contributions. 


8.5.1 Physical observables 

Let us consider the influence of these heavy top effects on the various LEP observables. 

Observe first that we could also have defined the p and t parameter from the Z —» p + p 
and Z —> bb width: reinserting the phase space factors, one has 


T(Z^p + p~) 
T (Z -»• bb) 




(8.67) 


fJtiV = 1 - 4s 2 , SM = 1 

4 2 

9bv = 1 _ 3 s + r 1 g bA = 1 + T 


( 8 . 68 ) 


3 In [106] an explicit check of the Ward-Takahashi identities used by us has been given 
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where 


s 


2 


1 

2 



Attol 


V2G ft M 2 zP/ 


(8.69) 


Formulas are greatly simplified by the possibility of writing the heavy top corrections as 
a redefinition of the Weinberg angle, through the relation in Eq. 8.69, and of the axial and 
vector couplings of the Z field with the b field. The results can then be included in the so 
called “improved Born” approximation for the LEP observables, together with part of the 
QCD and “pure QED” effects. 

For instance, considering the effect on the e + e _ —> /i + /i~ forward backward asymmetry 


AU - 

a fb ~ 


crp — ctb 

&F + &B 


where 


a F(B) 



(8.70) 


(8.71) 


in this limit it can be simply written as 


au — 

A FB — 


(?lv + gjLl) 


(8.72) 


and depends on p through the effective Weinberg angle s. 

But to illustrate the effect of this two-loop computation on the electroweak precision 
observables currently measured at LEP, it is necessary to include also the other relevant 
corrections: in Figs. 8.12,8.13,8.14 we report the Standard Model predictions for the Z leptonic 
width, r (Z —> e + e _ ), the forward-backward e + e _ —► p + p~ asymmetry at the Z pole, A BB 
(both affected by p only) and the total Z width, Tz (affected by p and r). 

The graphs are obtained using the ZFITTER code [71]; the dashed lines include the mf 
effects irreducible, that is, coming from the two-loop computations, for different values of the 
Higgs mass. 

Superimposed is the band of experimental values: 


-pTot 

= 2497.4 ± 3.8MeV 


pZ 

1 z 

= 83.96 ±0.18MeV 

(8.73) 

A ^ 

A FB 

= 0.0170 ± 0.0016 . 



To judge the importance of the two-loop effect, the dot-dashed curve at mn = 2TeV 
can be compared with the one-loop result, at the same value of the Higgs mass, given by 
the full line and resulting by setting to zero the second-order coefficients p ^ and r®. The 
rriH = 2TeV has been chosen because the effect is larger, but unless the top mass is pushed 
at unreasonably large values, the effect is at most of per mille. 
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8.5.2 Validity of the gaugeless approximation 

Is the gaugeless approximation sensible, for the real world? We have seen that for rrif in the 
150 — 200GeV range, the corresponding range of gt — 5.74 x 10 -3 mt/GeV is 0.86 — 1.15, 
while the gauge couplings are g = e/sin 6w — 0.64, g' = e/cos 9w — 0.34; alternatively, 
recall that in the two-loop computations we have systematically neglected M\y , Mz when 
compared with mt , a procedure which cannot be justified on numerical grounds. 

It is useful to consider one of the original computations at one loop of the heavy top 
corrections to the Z —> bb, performed in [41] without the M\y <C mt approximation, in the 
renormalizable t’Hooft-Feynman gauge. It is possible to separate the contribution of the 
graphs with the exchange of IV^ and the (fr^ fields and the result is that even for mt > 21% 
the “leading” ml contribution is only a factor of two greater than the subleading one. 

This means that while the large limit is sufficient for an order of magnitude estimate, 
the actual accuracy in LEP experiments requires the inclusion of the one-loop subleading 
contributions. 

At the two loop level no complete computation exists, but other authors [107] have at¬ 
tempted to estimate the subleading contributions by working in a reduced SM model, limiting 
themselves to an SU( 2) gauge group 4 . They find that the leading O and the sub¬ 
leading O corrections are of the same order of magnitude, and of the same sign. 

Similar results have been obtained in [112], estimating in a particular process (neutrino scat¬ 
tering on a leptonic target) the subleading terms, which are again found to be even larger than 
the “leading” ones. This results are a clear indication that at the moment the two-loop heavy 
top effects can be considered as a partial estimation of the theoretical error of the one-loop 
approximation, and that whenever the two loop heavy-top effects become observable, they 
will require a much greater theoretical effort to be compared with experiments. 


4 in other words, QED effects are not present 
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m H / m i 



m H / m t 


Figure 8.10: Asymptotic expansions and exact results for p^ and r® 
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Chapter 9 


Conclusions 


and perspectives 


We have presented in this thesis the computation of some radiative corrections to weak decays 
in the context of the Standard Model. We believe that the processes we have studied, the 
6 —> sy and the Z —► bb decays, shall receive further attention in the future, thanks to their 
sensitivity to New Physics. 

For what concerns the b —> sj decay, we have seen that the prediction of the inclusive 
decay is in good agreement with the experiment, and that this may be regarded not only as a 
success of the Standard Model but also of perturbative QCD itself: the dramatic enhancement 
of the rate resulting from the resummation of QCD logarithms is something which is difficult 
to miss. 

The study of the O (rrif) corrections to the Z —> bb decay did not result, on the contrary, 
in a dramatic effect. We have seen that the corrections to the one loop expressions are 
relatively small, improving our confidence in the phenomenological analysis of the LEP data. 
Our opinion is that these corrections are not sufficient to account for the bulk of the two loop 
effects: presently they can be regarded as an estimate of the error implied by the perturbative 
expansion. 

To conclude the work it seems interesting to us to give an example of a phenomenological 
analysis based on these two independent processes. 


9.1 An example: bounds on two Higgs doublets models 

We have seen in Chapter 2 how the presence of an extended Higgs sector may affect the 
prediction for the b —> s'y decay rate, in dependence of the mass of the extra charged Higgs 
particle, M H ±. It can be easily grasped, on the basis of the discussion in Chapter 6, that the 
same extended Higgs sector contributes also to the Z —> bb process. 

We have already discussed how this contributions are parameterized by two other quanti¬ 
ties, M h ± and tan (3: a comparison between theory and experiments allows to exclude portions 
of these parameter space. 

The limits coming from the b s'y process can be defined observing that in Model II the 
contribution from the extra doublet enhances the rate, for smaller values of M#±. We define 
excluded values of M H ±, tan (3 the ones for which the inequality 


BR th (B -> Xs'y) - 2 x 


theo. 

error 


> BR exp (B -*■ Xs'y) + 2 x 


exp. 

error 


(9.1) 
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Figure 9.1: Lower limits for the mass of the charged Higgs boson 


holds. This relation gives a lower limit on M#±, displayed by solid lines in Fig. 9.1, varying 
nit around the central CDF value. 

We refer to [109] (and references therein) for detailed formulas on the bound from the 
Z —> bb process: in the same figure the dashed lines display the resulting bound for the 
central CDF nit mass. 

For completeness the lower bound coming from the b —► ctv t [42, 98] is displayed with a 
dotted line. It is obtained by the simple formula 

tan/3 < . (9.2) 

CeV 

We can appreciate that in the analysis of this extended model the two processes are 
complementary, and they already exclude large portions of the plane. As larger values of 
tan /3 are preferred, in order to explain the large splitting in the (f, b ) multiplet, one can see 
that the b —► 57 process gives the most important constraint. 

9.2 Future work 

In the study of the b —* sj process, presently theoretical and experimental accuracies are 
comparable. We can expect in the near future an improvement of the statistics of the CLEO 
experiment, and therefore the b —> sy decay shall require a complete next-to-leading compu¬ 
tation. We have seen that such a computation is mandatory if one is willing to obtain sensible 
answers to the question whether the effect comes in its entirety from the Standard Model or 
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not. The techniques we have applied to this computation are in our opinion well suited to 
simplify the required three loop computations. 

Once a complete three loop evaluation of the QCD corrections will be available, it will 
be necessary also to re-consider the determination of the Effective Hamiltonian, obtained by 
the matching procedure at the higher evolution scale, as described in Chapter 3. In order to 
exploit the higher precision in the evolution coefficients, it will be necessary to extend up to 
two loops the matching, not only for the Standard Model but for all the extensions one is 
interested to compare, including the Minimal Supersymmetric Standard Model. 

For what concerns the Z —» bb process, and more generally the corrections to the observ¬ 
ables measured at LEP, we have mentioned some works where different parts of the two loop 
effects have been computed. A complete evaluation of the electroweak corrections does not 
exist, and we cannot say, at the present level of experimental accuracy, that such an effort is 
mandatory. 

We can just say that a complete analysis of LEP observables at the “per mille” level of 
accuracy shall require the full inclusion of two loop electroweak corrections. We hope that 
the techniques discussed in this thesis may result of help in this task. 
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Appendix A 

Notations 


In this thesis both the pseudo-Euclidean (E) notation and the Minkowsky (M) notations have 
been used, so it is useful to give here a translation table. 


A.0.1 Metrics 

The metric tensor is given by 


Minkowsky: a ■ b = rf iV a ^by 

Euclidean: a ■ b = 8^ u a^b v 

(A.l) 

where rf v = (—1, +1, +1, +1) so that 

a • b = rj IJV afjjy = a ■ b — a^b^. The transcription in 

the PE notation is given in the following table 


rpMO A rpE 

JU - iJu ^ 

x Mi _ X E 

df 1 = id ,f 
d^ = df 

(A.2) 

For a vector field in E notation the fourth component is purely imaginary 


A mo = -iAf 

A m = Af 


AM _ • aE 

Af 1 = Af 

(A.3) 

^0 i ~ <Lt 4 i 

s~iM _ r^E 

^ij ~ ^ij 


hence for instance 




g m g m ^ = g e^ (A . 4) 


Fourier Transform 

We adopt the following conventions for Fourier transforms of fields 

m = / < A - 5 ) 

as an example 

u{v-)lnv{p + )A^{k) 


(A.6) 
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with k = p+ +p~ incoming in the vertex is the amplitude for pair creation in electrodynamics, 
written in momentum space. Using this convention the expression for the four momentum is 

d 


P M = (E,p) = 


&' = 


0 

dx, 


*a-- iv 


d 


= |-»' v 


hence one has the correspondence 
jP = d» = ip** 


(A.7) 


(A.8) 


Spinors and Gamma Matrices 


The correspondence for spinors in the two notations is as follows 

= ij] E p) M = = p> E = iplA 

where gamma matrices, defined by 

lp.lv + Ivlp = 2 Ppv 


15 = -*7°7 1 7 2 7 3 = +*7o7i7273 = -^ppupal^YYl' 7 
1 obey the following correspondences 


1 M 0 = — iif 

7o M = +*7f 

15 = lf = 71727374 


l m = if 
if 1 =lf 


For the computation of squared amplitudes it will be useful also 
l\ =H ll = -70 = 7° 75 = 75 
707^70 = Ip 

while for chiral theories the expressions of left and right projectors are 
p l/r = 2 (1 + 75) 

so that one has two component Weyl spinors 

1 5 R = R Ri 5 = -R 

15 L = -L L75 = +L 

obtained from a 4-component spinor via projection 

R = Pr 4'> L = P L p 

note that one has 


= LR + RL 4>ip4) = L7 m L + Ri^R . 

Unpolarized density matrices, and Dirac equations are given by 

E spin U Q (P)up(P) = (~ l P + m ) (iP + m ) U ( P ) = 0 

E spin v a (P)vp(P) = {-ip - m) {ip - m) v{P) = 0 . 

+0123 = -1 = -2 (? 7 “- VpVp 


(A.9) 

(A.10) 

(A.ll) 

(A.12) 

(A.13) 
(A-14) 

(A.15) 

(A.16) 

(A.17) 

(A.18) 

(A.19) 
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A.0.2 Lagrangian and Generating functional 

The measure in the Generating Functional for Green functions is different in the two metrics: 
the correspondence is 



C E = _ C M (A.20) 


The generating functional for Green functions is given by 

Z [J] E = exp —S E + / dxJ ■ 4> Z [J] A1 = exp iS M + i f dxJ ■ 4> 
The generating functional of connected functions is given by 
W[J] E = In Z[J] E W[J] M = i In Z[J] E . 

Defining the “classical field” 4> c functional of J, as 




one obtains by Legendre transform the generators of 1 PI Green functions 


T [< j> c \ E = —W [J] E + f dxJ ■ 4> c T [( j> c \ M = W [J] M - JdxJ • <f> c 


(A.21) 

(A.22) 
(A.23) 
(A.24) 


A.0.3 QED in Minkowsky notation 

The lagrangian is given by 

Zqed = ~\ + m) ^ (A.25) 

where 

F IIV = d^A u - d v A il V = dfj, - i e Q A^ (A.26) 

note that the electric charge Q of the electron is —1. 


Charge conjugation 

The Dirac equation is given by 

[(dfj, - ieQA M ) + m] ip = 0 (A.27) 

taking the adjoint and transposing one has, equivalently 

—7/f + ieQA M ) + m 'ip' 1 = 0 (A.28) 

the charge conjugation operator C, equipped with the property 

C~ x ilC=- lti (A.29) 

allows to define the charge conjugated spinor ip c = C' _1 '0 T obeying the Dirac equation with 
the reversed sign for the charge 

[(dft + ieQA M ) ^ + m\ ip c = 0 . 


(A.30) 
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It is possible to choose C with the following properties: 

c T = c ] = -c = c- 1 

CC ] = C ] C= 1 

C 2 = -1 . (A.31) 

Some useful relations are 


c- l ijc 

= 75 


c^c- 1 

II 


crfc- 1 

II 

1 


crfc- 1 

II 

(A.32) 


A.0.4 SU(N) gauge theory in Minkowsky notation 

The pure gauge lagrangian is given by 

= < A - 33 ) 

where 


G% = d,Al d v A° + gf^^At 

and a = 1, • • •, N' 2 — 1. labels the generators, while the fermionic part is given by: 


(A.34) 


£ fermion — ^ U/?. ibj 1/)/) 

i=l 

where V M e is the covariant derivative. 

V tl = d tl - igT a A ( ^ 
and matrices T a satisfy the algebra: 

fj-ia r p^\ _ j jabcrpc 


(A.35) 


(A.36) 


(A.37) 


The field tensor G^ u can be expressed in terms of the commutator of covariant derivatives 
PV = -ig T aG% (A.38) 

fixing so the sign of the charge in A.34 on the basis of A.36. 


gauge transformations 

Some basical properties are 

- gfabc^A^ 6^ = -gf^co b G^ 
Sojip = iguj a T a il) Sui? = —igu> a il>T a 


(A.39) 
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Generators of the algebra 

Fundamental representation: 

N 2 - 1 

T a T a = C F I where C F = m (A.40) 

Adjoint representation: 

T b a c = ifabc (A.41) 

T a T a = C A I where C A = N (A.42) 

Minkowsky-Euclidean correspondence 

iS M = i j ( dx) M G tlv G IJV - $ ( 7 ^ + m) ^ = -S E (A.43) 

where 

S E = j ( dx) E Q G^Gp, + ^ ( 7 ^ + m) ^ (A.44) 
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The projection on the canonical basis of any product of gamma matrices can be obtained by 
using recursive relations 


-£(-!)• 

i 

E(-ir 

i 


(B.7) 

(B.8) 


. The following trace identity holds: 

Tr ( 7 M 7 W) = n\6 m ’ n s m g™ TV(1) 

n(n — 1) 

s n = (-1)—5— • (B.9) 

Some contraction identities: 


7ww«-A*n7A*i = (-!)” 1 (d - n + 1) 7 ^...^ 

7MlM2---Mn7MlM2...Mn = 

n—1 

Cn = (d- i) ; (B.10) 

i =0 

the c n coefficient can been expressed as a series in e defining 
c n = c n (0) + e c n (l) + £ 2 c n (2) + 0(e 3 ) . 

They are also useful the following commutation identities: 


{7/u ,M2v^n ’ 

hm 


lv} 


2 7/U,£i2,-..,Mn,Z' 

n even 

2 7/U,£i2,-..,Mn,Z' 

n odd 


(B.ll) 


and the reduction identities: 

7ju7 (n) 7/i = (-l) n (d-2n) 7 W 

n—1 

7 (n) 7/^7 (n) = (2n - <0 II (* “ d ) 7/x n > 2 


= (-1)’ 


2=1 

d — 2n 


d 


(B.12) 


1 


^W* ) - d(d _ 1) L 

7 W 7 W 7 W = / mi „ 7 <»> 

min(m,n) 

fm,n — ?Tl!5 n 5 m _|- n ^ ^ ( 1) 

1=0 


(d 2 -(d- 2n) 2 J 5^ + ((d - 2n) 2 - d) 7/i7 ^ 


7 W 7 ( n ) 


i I n 
l 


d — n 
m — l 


(B.13) 

(B.14) 


In Eq. B.12 summation over repeated indices is understood. 
Again the / coefficients can be expanded in series of e 


fm,n — /m,n(0) + £ fm , n ( 1) + 0 ( e 2 ) , 
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imposes the condition 

(Tr(l)) 2 = f^ j=2 rf . (B.18) 

It follows that Fierz transformations (B.16) are inconsistent with the usual choice Tr (1) = 4 ! 

On the other hand, as the values of f m .n follow from purely combinatorial relations in 
Eq. B.14, contraction identities remain valid. 

We have used in the computation reduction identities of the following form 

<8> = 7 < - n+1 ^ <8> 7^ n+1 ^ + n (d — n + 1) 7^ n_1 ' ) <8> 7^ 1 ^ (B.19) 

7 (n '*7 / u <8> 7^7/i = 7 ( - n+1 ' ) <8> 7^ n+1 ^ + n (d — ra + 1) 7^ n_1 ' ) <8> 7^ 1 ^ (B.20) 

7 M 7 (n) <8> 7 (n) 7 M = (—l) n ( 7 (n+1) ® 7 (n+1) - n {d - n + l) ^ n ~ 1 ^ ® (B.21) 

The relations in Eq. B.19 are fundamental and allow to reduce recursively any product of the 
form 

7w • • ■ 7wc7 (n) 7wc + i • ■ ■ 7 nm ® 7/mu ■ • • 7^ (h) 7 (n) 7^ (h+1) • • • 7^ (m) (B.22) 

where n G S rn is a permutation of indexes. 




Appendix C 


Standard Model 


In this appendix the notations for the Standard Model are set, in Minkowskian space, using 
Tj llL/ notation (see App. A). 


C.l Generalities 

The electro-weak model is based on a SU (2 )i e ft x U( l)y symmetry: the conserved hypercharge 
Y is associated with the 17(1) sector. 

C.1.1 Matter fields 

The matter fields are organized in SU(2) doublets with left chirality, and right singlets: we 
resume in the following the quantum numbers and the composition of doublets. 


Leptons 

L i = h ( 1 + 75) ^ j 


R l = \ (! -75)e l G {e, //, r} 


Assignments of weak hypercharge: 


Q = h + -/ 


YU = -UYRi 


-2 Ri 


Quarks 

where 

q\ € u, c, t qi € d, s, b . 

Weak hypercharge assignments: 

y Lq-f = g Tq T 1 


J R q r = 5 C 1 “ 75) 9 r 
1 = 3 (! “ 75 ) 




94 . ” 



(C.l) 

(C.2) 

(C.3) 

(C.4) 

(C.5) 
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The mass eigenstates do not coincide with the states that diagonalize the weak interaction 
matrix: the relation is given by the well known Cabibbo-Kobayashi-Maskawa matrix 


L' u = \{ 1 + 76) J 


d' — V ud d + V us s + V u bb 


L 'c ~ \ (! + 75 ) 


s' — V c dd + + V c bb , 


L t = 2 (! + 75) y 


b' — Vtdd + V ts s + Vtbb 


where the primed fields are eigenstates of W.I., while the non primed are mass eigenstates. 
The notation will be used for the elements of V, and V*^ for the ones of V'. 

C.l.2 Gauge field content 

Four gauge fields are associated to the local version of the SU (2) xU( 1) symmetry: three 2 3 
associated with SU{2) rotations, one £> /( with the remaining U{ 1) (hypercharge) symmetry. 
The respective charges are g and g’ , and the covariant derivative is given by 

V,, = Id tl - igl% ■ T - ig'B^ (C.7) 

while the classical lagrangian is 

71 SU(2)xU (1) — Bgauge + Bleptons T Bquarks T 71 Higgs T 71 Yukawa (C.8) 

and the Higgs and Yukawa terms will provide the mean to introduce masses in an invariant 
manner. 

The pure gauge lagrangian is then 


7 Igauge = ~ \ ^ 


where 


W IIV = d»W„ - cXAVg + gW„ x W v 


with (w„ x W v ). = eijkWiW* 


Bgy = d^By - dyB qi 


(C.10) 


(C.ll) 


C.l.3 Matter fields lagrangian 

The lepton lagrangian is 


7 '-lepton = - X - X ^7 • BRi 

l=e,n,T l=e,/i,T 

while hadrons are incorporated in 

71 quark = — ~ 5Z -^?t7 ' — 

q-^=u,c,t q^=u,c,t q^=d,s,b 


(C.12) 


(C.13) 
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C.l.4 Higgs fields 

The Higgs field is defined as a doublet of complex fields transforming under the fundamental 
representation of SU(2 ) with assigned hypercharge +1. 


4> = 


Y4* = 4 


(C.14) 


The form of the Higgs lagrangian is fixed by the gauge symmetry, while the couplings in 
the potential part are chosen in order to allow a spontaneous symmetry breaking. 


Cmggs = -V&V* - ^4* - A (V<h) 2 . 

The interaction of the Higgs field with matter fields is then of Yukawa type: 


(C.15) 


GYukawa — {-Yukawa—quarks T ^Yukawa— leptons 


(C.16) 


GiLi ■ <&Ri + h.c. 

Gf qi L’ qr -$R qi + h.c. 

G^ qi L’ q ^$R qi +h.c. (C.17) 

where the conjugated doublet 4> is given by: 

4> = ^ ^ Y <I = -4> (C.18) 


^-•Yukawa— leptons 
{-•Yukawa—quarks 


and the Yukawa couplings are proportional to the ratio of the quark mass to the massive 
gauge field mass 


Gi = 

Gf = 


Y2Mw 

9 

vT M n 


m 


m 3 V; 




— J 
T: V2M W 


m t 5. 




(C.19) 


C.l.5 Symmetry breaking 


Choosing a negative sign for the mass term, /j? < 0, the Higgs field acquires an expectation 
value: we may choose the phases in order to have it in the form: 


*' ( 75 <* + *) ) + ( £* ) 


where 

v 



(C.20) 

(C.21) 


(C.22) 
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Thus we obtain for the lepton lagrangian 


H)ll + i- 


£ Yu k awa= ~ "H {l + — H )U ^ — 

- —-rJ-—mi ( Ui(j) + Rl + T(jrLvi) 

y2 mw 

Analogously for the mass lagrangian of the quarks and their coupling to H, <1>( + ’ x- 

^-■Yukawa—quarks — £"mass,H T T , 

H - quarks lagrangian: 

£mass,H = ~ ( 1 + „ 9 H) V m q qq ; 

V 2 m w J 


(C.23) 


(C.24) 


(C.25) 


quarks 


■* - quarks lagrangian: 


£*= - A^(R^ + 4 +g' 4 r«i) 


y/2 mw ' ' 
g ml. t j 

+ — j= -( Vj 

\f2m w V J 

X - quarks lagrangian: 


X0 U R + Vij ul R ( t >+ d'L 


Cx = i 2M^ x ^ (^9175^-^^75^) ; 


(C.26) 


(C.27) 


quarks 


Weinberg angle 

The gauge covariant interaction of Higgs and gauge fields gives rise to mass terms for the 
latter: one gauge field remains massless, corresponding to the unbroken U{ 1) electromagnetic 
symmetry. To put in evidence this symmetry breaking pattern it is convenient to define 
“rotated” fields 


W'l = Z jL cos 0\y + An sin Ow 
B M = A tJ cos 9 W - Z M sin 9 W 


(C.28) 


A^ = B fi cos 9w + H 7 ^ sin 0 W 
Z lt = W^ cos 9w — Bn sin 9\y 
The covariant derivative turns into the following 
= Id,, - ig ( W^T 1 + Wf L T 2 ) 

- i cos 9w + gT 3 sin 9 W j A M 

- i ^ gT 3 cos 9 W ~ g'\ sin 9 W ^ Z M 

The condition that the field A /t is associated with the photon imposes that 
g'\ cos 9w + gT 3 sin 9w = e (t 3 + = eQ 


(C.29) 


(C.30) 


(C.31) 
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e = g cos 6\v = g sin 9w 


(C.32) 


hence the relations 

j 


g_ 

g 


= tan 9w 


e = 


g'g 


vV 2 + 


g- 


(C.33) 


We may then define charge lowering and raising operators of the weak isospin (and hence of 
the electric charge, according to Eq. C.31 


T ± =T 1 ± iT 2 . 

We may individuate the charged components of the gauge fields: 
1 


(wl^iwl) 

and the covariant derivative can be simply written as: 


W = 

11 V2 


^ w l 9 

v, = w„-- 2 


( W+T+ + W~T~ ) - ieQA^ - zeQ'Z,, 


where the neutral charge matrix Q' is introduced 

Q' = T 3 cot 9\v —— tan 9w = —— j, -— 

2 sm 9]y cos 9\v 


(T 3 - sin 9y V Q 


(C.34) 


(C.35) 


(C.36) 


(C.37) 


C.l.6 Field transformations 

It is useful to reformulate transformation laws in terms of rotated fields: to this end let us set 


J± = 72 ^ ^ iuj2 "> 
u A = 7 cos 9\y + w 3 sin 9\y ljz = w 3 cos 9\y — oJ sin 9\y 

It follows for the transformation of gauge fields 

5A^ = d^ A -ie(uj~W+ -oj + W~) 

SZfj, = d^uz - ig cos 9 W (u~W+ - oj + W~ ) 

5W± = =f ie (w ± d At - lu a W =f ig cos 9 W - oj z W^ 

while for the Higgs it results 

5H = ~ (oj—(j)+ — u>+(f>-) + ojzX ■ L 
2 sm 2 0w 

g 6 

= -(uJ-4>++uJ + (j)-)- ujz(H+ v) . 

Z Sill Zu\y 


(C.38) 


(C.39) 


5(j)± = ± 


l { u± {h+v) - 7 ±x ± ie ( wa + 


and finally for any doublet or singlet of matter fields 


<54' = 


— (uj + T + + u> T ) + ieQojA + ieQ'uJz 

v2 




(C.40) 


(C.41) 
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C.l.7 The complete Lagrangian 

The gauge lagrangian may be expressed as follows 

C nmme = --W±W~^ - -Wlw^ u - -B,, V B^ U 


where 


wf UJ = d»w*-dvwl-ig(w+w;-w+w-) 

that is, evidencing the subgroup 17(1) of electromagnetism 


(C.42) 


(C.43) 


Bgauge ^ B U1 J F ahi ^ZgyZgy 


+ igfawf-dyWfiw+W- 
+ 1/ (w+w- - w-w+) w+w- 

- \\{d,- igWl) W+ - (r% - igW't) W+ 


where 


Bfxp — d^Ai/ d v A^ 


— dgZ v — dyZg, . 


The hadronic lagrangian is, in full form 


f~-quark — ^ i&QA^ ieQ Z q 

quarks 

+ -% (Vji<'fw+A + vp L rw- .u{) 


(C.44) 


(C.45) 


(C.46) 


Separating the Higgs lagrangian in kinetic and potential parts, CHiggs = 7 Higgs ~ VHiggs we 
have 

^Higgs = ~ + + i-^jzzW^ (f>* ^ 


x (Vc(> + -i-ir 3 ^ + -i-B^$ + - i-^W +ti (j) 0 ) 

V 2 2 \/2 / 

x (V</> 0 + *|lT 3/ Vo - - i^=W-^ + ) 


(C.47) 


that is, putting in evidence the charged gauge fields 


7f uggs = ~ [d^ + i g -W*4> + i^Bn4> 


x >+ - i|fT 3 G/>+ - i^B 11 ^ j 




C.l Generalities 


122 


x (d^o + i| W^cj> 0 - i S -B ^ 0 ) 

+ i-^=W +fi (</> 0 <9 m (/> - - (t>~d^(t )o + ig'B,^ 0 <j)~) + h.c. 

_ X -g 2 W^W~ ( 0 + 0 - + <f>* 0 <f> 0 ) • 

It is convenient to define the covariant derivative in the subgroup U( 1) x 17(1): 


/’3rri3 ■ ^ 


B fl = Id, - igW'ijT J - ig — B fl 
= dp — ieQA, - ieQ'Z, 

Q = — -( T 3 ~~ sin^/ Q) 

sm^cosw v ' 

An useful relation is the following 


2 

Setting 


(gWl- jB^ = Y—Z, = m z Z, 


v_H + ix v 

V2~ V2 v/2 


V„s- (v„sy =iV2^d F x + 


- 3 Zz„H 

2 e M 


one is able to rewrite the kinetic Higgs lagrangian in the form 

Tm ggs = ~ (v^+yv^ + - (V,S)*V^S 

- razZp^x+l^-Z^H) 

+ imy/W^ (<9 M + ig'B^ ) <f>~ + h.c. 

+ Aw+ (Sd' l 4>~ - (jrd^S + ig B' L S(f>~) + h.c. 
v 2 

- 7 2 9 2 W+W-» (</»+</»- + S*S) - gmwW+W-"H 

- rn 2 w W+W-' 1 - ^rn 2 z Z lt Z' 1 . 

Noting that 

\ — l r ,2 ( mfl ) 2 /; 2 __l rr) 2 

A - 89 \mw) V - 

it is possible to rewrite the potential part as follows 



(C.48) 

(C.49) 

(C.50) 

(C.51) 

(C.52) 


(C.53) 

(C.54) 

(C.55) 
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C.2 Quantization 


A perturbative evaluation of matrix elements requires to break the gauge invariance and 
introduce a gauge fixing. In this section we want to discuss a few technical points that result 
of great practical advantage. 

C.2.1 Non linear gauge fixing 

After the symmetry breaking, the interaction between gauge fields and non-physical Higgs 
fields gives rise to spurious terms in the quadratic part of the lagrangian: it is possible a gauge 
choice which explicitly eliminates this coupling, avoiding so mixed propagators between gauge 
and Higgs fields. The gauge-fixing is given by: 

Cg.f. = ~ {d^f - i- (d,z^ + gm, zX ? - j |(^ - igW^W^ + i^m w <t>+ | 2 (C.56) 

Putting together gauge fixing and Higgs lagrangian one has 


+ 


(AA)* - (v^s)* TPS - |m | x 2 - 

1 g£. 

2 e 

gg' 


-mzZ^H 
-m w W+Z^(j)- + h.c.) 


+ i-j=W+ (Sd^(p~ - (jrd^S + ig'B^Scj)-) + h.c. 
V ^ 

- \g 2 W+W-» (</)+</>- + S*S) - gm w W+W-' J H 

2 


1 


- rniyW+W-i- 1 - -rn z ZijZ 11 


At ■■ 2' 

We may collect terms in the total lagrangian containing only gauge fields 

£•gauge total 


(0.57) 


\W" - A 

11 1 

— Z^Z^ V — r7 ^ 2 ™ 2 


2g 


(dgZ^Y - -m 2 z Z,Z» 


+ 

+ 


ig (dgW^ - dvWf) W +fi W~ u 

\y (W+W- - wywy) w + »w~ 


(A - igW3) W+ - (d v - igWl) W+ 
A - igWfiW+tf - m 2 w wyw~ » 


/z / r r | W r T [i 

and terms coupling gauge and Higgs fields 

gauge—Higgs = ~ (A</> + ) ^A ~ (A* 5 ) ^S - |m|x 2 - 


(C.58) 
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- ^mzZ^H - gmwW+W-^H 

+ (^rnwW+Z* 1 ^ + h.c.) 

+ i-^W+ (Sd^cj)- - ct>~d^S + ig'B^S4>~) + h.c. 
v 2 

- + S*S) (C.59) 

of course the unitarity is reinforced introducing anti-commuting scalar fields: the ghosts. 


C.2.2 Ghosts and Slavnov Invariance 

It is convenient to rewrite the gauge fixing lagrangian in the form 


■g-f- 

= ~—F\ - —F| - -F+F_ 

2a A 2r/ z £ + 


Fa 

= 


F z 

= dgZ 1 '- + gmzx 


F± 

= (d,, hF igW 3 ) ± itmwQ* 

(cm) 


The properties of transformation under SU{2) x U( 1) are 


5uF A 

5uF z 


S*F± 


nu A - ied^ (u~W+ - w + W“) 

Ou>z — ig cos 9\vd^ 1 — u + W~^ 


rjrriz 


ojz (H + v) 


— C — -UlU-$ + +LU + &-) 
sm 2t) w 2 


d» (d lt T igWfj =f igWl (d> 1 ^ igW 3 ) 

iguj 3 (& 1 T igW 3m ) W± ff g 2 W± (u~W^ + - u + W^~) 

uz ) ^ >± 


- £rn w 


l“ ± (HFv±ix) + e(u A + tan2ew 


1 


(C.61) 


Following Becchi-Rouet-Stora we may define the gauge transformation with parameters w* = 
X l (, iJ = X'C, where the compensating helds X®, X' anticommute, so as the x —independent 
parameter Q. The BRS transform s is defined as the right derivative with respect to £ of a 
gauge transformation parameterized in terms of fields X. In the basis defined by the same 
rotation which diagonalizes the mass matrix of gauge fields one has 


X ± = j, ( x 1 T iX 2 ) 

X A = X' cos Ow + X 3 sin 9w Xz = X 3 cos Ow ~ X' sin 9w 


(C.62) 


and this fields are introduced in the action with the following lagrangian 
Lg. c . = + X a sFa + XzsFz + X + sF + + X sF L 


X a sFa= - & J X A d lx X A + ie&‘X A (x-\V+ - X+W~) 

X z sF z = - fFXzdgXz - vm 2 z X z X z + ig cos d w &*X z (x~W+ - X+W~) 
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X ± sF ± = 


- r,m.z^XzX z H + V m z | (. X Z X~<!> + + X z X + <5>~) 

- (&> db igW^) X ± (d„ hr z<?W 3 ) 

db igX^X 3 (V ff i^W^ 3 ) B 7 ^ - g 2 (x^X^W* - X ± X ± W^ W^ l± 

- £m w 9 -X ± X ± (H db z X ) - tm w e (x±X A + ****) 

(C.63) 


where s is the BRS transformation. 

In the following for completeness we report the effect of BRS transformation. 

The explicit expression for gauge fields is 

sA» = d lx X A -ie(x~W+-X + W~) 

sZ „ = d ll X z -igcaad w (x-W+-X+W-) 

sW± = dpX^^ie^Ap-XAWfj^igcxxOw^Zp-XzWf) (C.64) 

while for Higgs fields 

sx = f (x-i»+ + x+$~) - X^x z (H + „) 

*** - ±2 X *(ff + ,)-|x^ ±fe (^ + i= ^ft)^ (C.65) 

and lastly for the compensating (ghost) fields 
sX i = i e ijk X j X k sX B = 0 

sP = |F ?: sA B = (C - 66 ^ 

that is, 

*** = Tie (v, + ^X z ) JC± = IF t 

sXi = -ieA+A" sA a = ±F a (C.67) 

sX z = \F Z • 

It is easy to verify, thanks to the nilpotence of the BRS transformation (s 2 = 0), that the 
variation of the Gauge Fixing Lagrangian is compensated by the variation of the Ghost terms. 
The original gauge invariant lagrangian is also BRS invariant, hence at a classical level the 
gauge fixing breaks the classical invariance but leaves a BRS invariance. In App. D we shall 
see how it is possible, and in many cases convenient, a different choice of the gauge fixing 
which enables to preserve also the classical invariance. 


C.2.3 Feynman rules 

It is now useful to report the Feynman rules for the electro-weak sector of the Standard 
Model, for the given choice of gauge fixing. The modifications required by the Background 
Field Gauge will be discussed in App. D. 

The convention chosen in the following is that momenta and charges enter the vertex. 
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Masses 


W ± (Z„) 


m w (m z ) 

\flm w 

y/rjmz 


Propagators 


scalar field 


fermionic field 


vector field Z 


vector field A 


(C. 68 ) 


(d* (-p) $ (p)) = 


p 2 + m 2 


{*M>) = ~ 


i/p + m 


vector field IT± (-p) W? (p)} = -(!-£) 


(Z„ (-p) Z F (p)) = j—— 7 [i)„„ - (1 - if) 7 

p 2 + m 2 |_ p 2 + i/m 2 J 

(-p) 4, (p)> = ^-^2 ^ -(!-«) p /^ A2 


Matter fields-gauge fields couplings 


(C.69) 


(uiW+i,) = 

(iw-uj) = -^= 7,7 

(p,HV+e) = 
i'ipqAfj.ip) = —eQq'y^ 
(uZ^u) = —:--- 


sinu/cosw l_2^ M 


■ sm w 7 ^ 


^ = ^ ein^eoSH.- [-pe+r in2 "' 7 >-. 

<(Z„1> = el [tan,,/ ( 27 * + 7 J) - cot,,- 7 J 

{uiZ^i) = [tan w + cot w ] 

Matter fields-Higgs fields coupling 


(C.70) 




*9 ,/ w 

3 I 

V 2 V mvK 
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(<k® Uj ) 

{m + l) 

(U~Ul) 

{fHf) 

(hi) 

(uxu) 

(dxd) 


V2 


V 






-ig m 

/o ^ ''' 

V 2 rnw 
-ig mi 
v/2 m-vu 
i_9_ mf 
2 mvu 
5 m 

-«- 75 

2 myr 

, 9 m u 

+ 77 - 75 

2 mvK 

9 rn d 

-«- 75 

2 mvK 


Gauge fields self interactions 


^ j W+(p+)W~ (p-)j = 

i( e 

\ gcos w / 

+ f 1 - 7 ) (W - 


(p+ - P-)a Vpv + 2 Vp° ip+ +P-)p- 2r W (p+ + P-) e 


(w'+w-w+w^) 

(ApA^W+WJ.) 

(z,z,,w+ + w u -_) 

(A^w+w-) 


ig \ 2 Vpp'Vw' VpvVf.i'v' vpu 1 vu\ 

1 (~9 2 ) [ 2 9pp'9v+v- ~ Vpv+Vp'v- ~ 9pu~Vp'u+] 
l ^ ® ) [27/-i/i'7^+!/ _ — Vpu+Vp’u~ Vpi/~Vp'u+\ 

i(-g 2 cos^j [2g^g v+v - - V^+Vp'v- ~Vpu-Vp'u+] 
i(-egcos w ) [ 2 Vpp'Vv+„- ~ Vpv+Vp'u- ~ 9pv~9p'v+] 


Gauge fields-Higgs interactions 


(C.71) 


(C.72) 


($(p+)A^(p_)) 

($(p + )Z^(p_)} 

(W+Z^~) 

(ZpZ^H) 

(W*W-,H) 
(U) (p")) 


ie (p_ -?+)„ 

- (c/cosyy -</ sin W /) (i?-P+) 
.gg' 

% ^ hflwh\nu 

i gg' 

~2~ mzV ^' 

-igmwPpp' 

±| (p?-pj) 
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( p T ) x C p x 


z, 


( P Z ) Hx 


ZfiZ v 


HH 


w+w v 


XX ) / 

s XX J / 
{A^A ^ + $-) 
(ai^,$+$>-) 

(Z^Z^*-) 


(w+w M -<&+<&-) 

(w±7U7<f>^) 

(w^z^) 


I fe - p£) 

^ (5 cos^y +c/ sin W /) P,f 

(ff cosy/ +£f' sinjy) 2 T) IIU 


.1 2 

-*,5 tyu/ 

—i2e 2 ?7 A j A j/ 

—ie 2 (coty/ —tanw) Vh/j. 1 
-ie 2 - (cotvu -tanw) 2 Vnn 1 


■ 9 

1 2 'W 

eg 


2 

^.eg 

9 




V»v 


Higgs fields self interactions 


r 


(HHH) 



(HHHH) 

(HH X x) 

(xxxx) 





myy 


-i-gm w r 


-i-gm w r 


.3 
2 ; 

.1 
2 ; 

■ 3 2 2 

-Z-G ?’ 

4 

. 1 2 2 
-z-5 r 

.3 2 2 

-z-5 r 

.1 2 
-i-g 
4 


-i-q 

T 


2 

2 


(C.73) 


(C.74) 


Ghost masses 
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X z ■ \frjrn z 

X ± : V£m w 

Ghost-Gauge interactions 

(Xa ( P A ) X*W?) 

{-Xz \p Z ) *=FW± ) 

(X±X 3 W± (p±)) 


Ghost-Higgs interactions 


(- X Z X Z H) = 
{Xz^X^) = 



{X ± ^ ± X A ) = 
(X^Xz) = 


= Xiep^ 

= =F ig cos w pZ 
= ±ig (v - P)^ 

= -ig 2 pfiv 

= Tigp± 


—irjmz -— 
e 

g 

irpriz - 


—i£,mwe 

—i^mwe 


1 

tan 2 9w 


(C.75) 


(C.76) 


(C.77) 




Appendix D 

The Background Field gauge 


The background formulation of field theory is an useful tool, which allows to take full advan¬ 
tage of local symmetries. The point is that one can separate the field in a “classical” (we will 
see below in which sense) and a quantum part; to be definite one finds that for gauge theories 
the original symmetry is implemented via two different transformations, and the “classical” 
symmetry can be implemented regardless from the quantization method. 

In the following we will summarize known results about this formulation and its imple¬ 
mentation in the Standard Model (see for instance [110] and references therein). 


D.l Generalities 


To be definite, let us consider the theory of a single scalar field (j>, described by some lagrangian 
C ((f), d^,(f>). The Green functions of the theory can be obtained by the generating functional 


Z [J] = e lW ^ = j D^exp i j dxL + J<f> 


where W is the generator of connected Green functions. Let us consider instead the same 
theory in presence of a fixed background field, that is, 


e iW[<l> c \j] _ J £)(^ ex p ,[ J dxC (tf) + (j) d j + J(f> 


(D.l) 


Suppose then to choose J as a function of 0 cl , in such a manner that the condition 


d 


-W 


id 


, J = —J(x) 


d(f> cl (x) 

is enforced. Explicitly one has then 
dW f , dJ(y ) dW 


d(f> cl ( x] 


f, oJ (y) dW / (C] \ 

+ l dy ww)dm = ~ ^ 


(D.2) 


The sense is clear in the zero loop ( h —> 0 ) limit: in this case the expression for W is the 
extremum of the integrand under the variation of (f >, that is 


W 


<f) d ,J = extr ^ J dxL (<j) + + J(x)(f>(x) 
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which means that a solution is given by 
<f> = 0 


5C 


Jo = — 


W 0 I 0 c \J U c1 


acI 


= J dx£ 

Now equation D.2 can be iterated 

J = Jq + hj\ + h 1 J 2 + •.. 


(D.3) 


w = w 0 + nw 1 + h 2 w 2 +... 


J\ 


d 

d(f> c 


W] 


<t> c \Jo +o(ny + ... 


(D.4) 

The important result is that the functional W | 0 cl , J so obtained concides with the 

generating functional T (V/> cl j of 1 PI Green functions of classical fields. To see this, we have 
the give a sense to the term “classical” we have been using until now: we will choose the 
condition on (j) cl 


J 


cl 


= 0 


which substitutes the classical motion equation 

SS n 

~ J 

with the perturbative form 


lC l 


Sc/)' 

By now it is almost evident the thesis: in fact 
• shifting the fields one has that 


(D.5) 


(D.6) 


W 


A cl 


, J =w[J\- / J4> 


cl 


which means that one has the Legendre transform of W [J], that is T 
• The equation of motion implies 

so the functional T is an effective lagrangian containing the quantum fluctuations, which 
allows to compute Green fucntions by a tree level expansion in terms of irreducible 
quantum kernels. 

What about gauge theories? One has the following relevant consequences: 
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1. the gauge fixing needed to compute quantum fluctuations is largely arbitrary: by choos¬ 
ing it invariant under gauge transformations of the classical fields one ensures that the 
effective lagrangian F is invariant under the same classical transformations. 

2. As the amplitudes are built of irreducible functions derived by T and connected by clas¬ 
sical lines, the gauge fixing needed to determine the classical propagator is completely 
arbitrary, in the sense that any gauge dependence in it is canceled by the effect of Ward 
Identities for the irreducible vertex parts. 

Having laid down the basics, we will work now on our preferred model. 


D.2 Compact Notation 

We will follow the notation used by Van Damme in [24, 31]; let us first consider the gauge+higgs 
sector: 


C = 


~ ( D ^) t (D^) - V (V<h) 


/oa 

U [1V 

D,G a u 

(T a ) t 

rpa rjnb 


abc/^b r^c 
\ J,V 


= - d v Gl + / 

= d^G a u + f ’ 

_ rj-ia 

_ j^abcrpc 

= d u &-G“&T a 

V [A 


abc/^ia s~ib 

[l^V 


(D.7) 


where the index a will run over the Lie generators of the symmetry group, even for non-simple 
groups as is the Standard Model case: all couplings are contained in the definitions for the 
generators itself. The field $ will in general be a multiplet of complex fields. Simmetry 
transformations are of the form 


5G; = -d,A a + f 
<54> = A a T a <Y 


abcpbQC 




where in our notation A are real parameters. 

We want now to perform the shift $ —> 4? + <h cl , G b 
to show then that 


(D.8) 

Gu + G y 1 . It is a matter of algebra 


G a 
[IV 

D*G a v 


D*G a v - D + f abc G\G c v + G a b cl 
d„Gl + f abc G b u d G c v 


' [IV 


- (£>„<&)+ (D^) 




-G a „ \ ( D £ + $ cl ) ) 1 T a + 4* 


cl 


+ 


- (4> + 4> ci y T a D c j($ + <Y 

^ T a T b + $ cl ) 


(D.9) 
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where we have put in evidence the quantum gauge fields. 

Suppose now that the potential V is of the Higgs form, thus giving rise to a symmetry 
breaking: in this notation we may then write that the classical field acquires an expectation 
value, thus giving rise to a shift 

<h cl = d* cl + e . 


This gives rise, thanks to the last term of D.9 to a mass matrix term for the gauge fields: 

M ab = (T a T b + T 6 T a ) £ . (D.10) 

The Goldstone theorem states that this mass matrix will have some zero eigenvalue, so at 
least one ( as in the Standard Model ) gauge field will remain massless: this will enforce us 
to introduce some gauge fixing to work in perturbation theory. 

They appear in the quadratic part of the lagrangian also a set of terms involving gauge 
and Higgs fields, which would lead to mixed propagators: 


-g; 


d^T a e - e^T a d^ 


(D.ll) 


The resulting complications in the Feynman rules can be avoided by choosing, following 
t’Hooft as we did in App C, an appropriate non-linear form for the gauge fixing term: 


Gg.f. = ~ 


1 

2 ? 


D^Gl + 






(D.12) 


as the classical field acquires a “vev” e, the resulting quadratic term combines with the one 
in D.ll to form a total derivative. It is important to observe that this gauge fixing is invariant 
under the so called type II transformations 


6G,, = A A G lt <5G d = A A Gf A - <9 /t A 
<5$ = A a T a $ (5<h el = A a T a T cl 


(D.13) 


this results in the classical invariance for Green functions of classical fields: in fact the ac¬ 
companying transformation of quantum variables amounts to a change of variable in the 
functional integral. 

Finally after the shift the mixing between non-physical Higgs and gauge fields disappears, 
and the following quadratic lagrangian results for the gauge fields 

£ 2g au g e = (d,Gl - d v G*) 2 - ^ ( d.G*) 2 - ±G;G b Ma b (D.14) 


while for the (quantum) Higgs fields one has 

£ 2 hi ggs = - (5 /i d>) t - l -i JVr a £ - £tr°$] 2 - h 2 (e, $) 

and the last term results from the symmetry breaking potential. 


(D.15) 
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D.3 Application to Standard Model 


Let us write down a transcription table to use the results of the preceding section in the 
Standard Model. 

The following correspondence holds: 


G*T* = —igWfj • | - ig'B^- 


(D.16) 


giving so f abc = ge abc for {a, b , c} £ {1,2, 3}, and zero otherwise. 

On the other hand it is easier to perforin already at this stage the Weinberg’s rotation 
which gives diagonal mass matrices, and the choice of charged states; so we use the following 
“vector” of group generators: 


{T a } = j--^=<x + , --^=cr , ~ieQ, —ieQ' 


(D.17) 


where cA are the raising/lowering operators in the SU( 2) group, and Q,Q' are defined as in 
Sec. C.1.5. 

We need then to write down the gauge fixing term in this notation: first note that the 
“vev” has the form 


so it results that 


T+£ = ^ J T~£ = 0 

(-ieQ) £ = 0 (-ieQ') e = 2sin ^ e COSw £ 

Next let us write down the type II covariant derivatives for the gauge fields: 

D$w± = d tl wt T ig (w™W± - wpp) 

= d^Wt =F ie (A*wt ~ A u W± d ) 

T igcos w (z c JW±-Z„W± cl ) 

D*A V = d li A v -ie(w+ d W--W~ d W+) 

Dpu = d,Z u -igco Sw (w+ cl W~-W~ cl W+) 

We will need also the gauge fixing term, which is given by 

£g. f . = + [ 

- 1 [Dpp + e (-ie) (V qq* _ 4>t cl g<h)] 2 

- ± [Dp, + £ (-ie) (Vq'<*> c1 - 4>t cl Q'<h)] 2 . 


(D.18) 
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To rewrite this expression in terms of component fields, we need the value of the following 
terms 

a -$ci _ $tcl a - $ = _L _ * x ) _ 0+ (if cl + v - 2 X cl j" 

_ $ tcl Q$ = 0+ cl </,- - 0" cl (/> + 

- T tcl QM> = ---- \(l - sin^) U-cP +d - 0- cl 0+) 

sin M /cosw L\2 / v 7 

+ |((f c1 + t) X -^ X c1 )] . (D.21) 

To write down the gauge fixing lagrangian in expanded form it is easier to consider sepa¬ 
rately the different terms 

Quadratic part 

as expected, we have 

£g.f.2 = -| \^W+ + ^ M w4> + 1 2 + \ (d^) 2 + i + ^Mzx) 2 (D.22) 

Trilinear in gauge fields 

£g. f .3g = [w™ - r%lT-lT+) - (^1T+1T- Cl - d v w-w^) 

-f (w") (w- d W+ - W+ d W") (D.23) 

Quadrilinear in gauge fields 

£g.f. 4g = -y “ ^^ Cl ) ( W » dw » ^W- d ) 

_ e 2 +g 2 c° Sw (yy+ciyy- _ w -d w +^ (w-« w + _ W^W") (D.24) 

Trilinear in Higgs field 

£g.f.3H = 

e| (y) [# (V +Cl <r + <r c V + ) + (V + <r cl - <r</> +cl ) 

+Cx-—— gV ) *X [(l - 2sin^) u~(p +d - ^ + 0 _cl ) + i (,H cl x ~ H X cl ) 

2 COSJ4/ V2cOSvy / L\ /V /V 

(D.25) 
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Quadrilinear in Higgs field 
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g.f.4H 
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1 2 


1 — 2 sin 


w 


> +cl - t + r cl ) + i (h c \ - H X cl )] 2 (D.26) 


D.3.1 Modifications to Feynman Rules 

By expanding this expression in terms of component fields we obtain, as already anticipated 

• a quadratic correction to the gauge lagrangian which cancels the Higgs-gauge mixing 

• a gauge dependent correction to the cubic and quartic coupling of gauge fields ( one or 
two classical legs ) 

• a gauge dependent correction to the mass matrix of the would-be Goldstone bosons 

4^,X 

• gauge independent corrections (proportional to the “vev” ) to the cubic interaction of 
gauge fields and Higgs fields 

• a gauge dependent correction to the interaction of Higgs fields, (one or two classical 
legs) 

In this Appendix we shall limit ourselves to the vertices relevant for the calculation Z —> bb. 
Then we need to reconsider only the vertices of the y cl with 2,3 quantum Higgs fields. 


Relevant x cl -quantum vertices 



(; x cl x4> + (t> ) 
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2 cos 


i 

--^gmw 

-2i\ 



gy \ 

cos w ) 
— m z 


cosw 
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It results that only the trilinear vertex H , y, y takes a gauge correction. For completeness we 
list in the following the other relevant vertices, which are the same for quantum or classical 
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fields 

(xtt) 
(xbb) 
(U + b) 
(b<t>~ t) 
{tHt) 


9t 9 rn t 

-^15 = A-75 

V2 2 mw 

9b r = 9 rn, b 

x/2^ 5 2 mw 


mPi = 


igtPR = 


wt 


ig m t 
\/2 mw 
ig m t 
y/2 mw 
ig m t 
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v/2 2 m w 


(idem for b ) 
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